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ABSTRACT 


A linearized mathematical model-is developed which simu- 
lates the dynamic response of the Navy F/A-18 for the study 
of flight control reconfiguration. JThe airerateeesemoqeuas 
as a multi-input multi-output, sampled data, closed system, 
which couples the dynamics of the flight control system to 
the aircraft linearized small perturbation e quale tones The 
discrete time, state variable equations for the system are 
then formulated. <A computer program is developed which will 
compose the model matricies and compute the response of the 
G@licraft to Stiek and piidder anpues: 

To study flight control reconfiguration, the Model alee 
individual actuation of either a left or right control sur- 
face. Aircraft response to the actuation loss of either the 
left or right stabilator is simulated in the progran ea 
program is designed to implement the reconfigurable control 
Mixer, Cunrently Under Sean for the Self-Repairing Digital 
Flight Control System. 

The computer simulation was written in VS FORD RAT 
copy of the program and simulation results are included 


the appendices. 
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Le INTRODUCTION 


Tactical aircraft face airborne and ground based threats 
which continue to grow in number and capability. in fame. 
combat engagements the lethality of these systems will pose 
a considerable threat to sea survivability. In addiiiaaam 
to combat losses, the loss of aircraft in battle damage re- 
pair, or awaiting repair, will significantly redttce our tac- 
tical forces. Projections on the survivability OFM 
forces [Ref. 1] during the initial days of engagement indi- 
cate 68% of the tactical aircraft will be out of action 
after the third day of battle. Of this total, 22/2735 
lost in combat and 464 lost while in repair or awartingagee 
pair. This is indicated in Fig. 1.1 which comes from Kefer= 
ence 1. It is clear from projections such as this, Tiaigue 
PeGuees sac to remain a superior threat to the enemy, COmmaae 
ued emphasis must be given to reduce the combat vulnerability, 
and increase the reliability and maintainability of owmeseenee 
ele alk a ame maine, 

To improve the combat effectiveness of future tactveam 
aircraft, the United States Air Force initiated the Self 
Repairing Flight Control System, Reliability and Maintain= 
ability Program. Reference 2 outlines the program plan and 
goals. An Air Force sponsored study [Ref. 1] showed that 


Significant improvements in aircraft survivability and 
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BO% 
68% 
AIRCRAFT 
OUT OF ACTION 
60% 
40% 
/ oe es IN BAT Ge 
7 DAMAGE REPAIR 
7 
7 Roe 
7 ate 
7 eee ee ee PRIVION 
20% oe : ‘ . mo +? 
13% AWAITING 


MAINTENANCE 


DAY } DAY 2 DAY 3 DAY 4 DAY 5 
ASSUMPTIONS 


e 288: ATTRITION | 
© 4701 DAMAGE LOSS RATIO 
e USAF DAMAGE REPAIR TIMES 


Figure1.1 Tactical Aircraft Battle Damage Repair Statistics 


i! 


reliability can be achieved with a self repairing digital 
Plivchtae Otero ts ys teme 

One of the principle features of a self repairing system 
is flight control reconfiguration. A reconfigurable fiom 
control system automatically counters aircraft loOSS Of See 
trol due to impairment, or loss, of a control effector . eee 
fundamental idea is to utilize the existing undamaged effec- 
tors to preserve the normal flying qualities of the unimpaired 
aircraft. The approach outlined in Reference 1 is to design 
a reconfigurable control mixer to be placed between the flight 
control laws and the control effectors. In the face of air- 
Crate impairment the mixer would reallocate control commands 
to the unimpaired effectors so that flight critical pie 
yaw, and rolling moments would be preserved to the maximum 
extent possible. Using the control mixer concept, the exmeee 
ing flight control laws would not be altered. 

Although the development of the self repairing system is 
intended for the advanced tactical aircraft, it is peOSgagimae 
that the mixer could be implemented in existing airframes. 
Tae could be done with a control reconfiguration module 
interposed between the aircraft flight control computer cea 
fly-by-wire actuators (her. Jae 

For this thesis the McDonnell F/A-18 was chosen to study 
the reconfigurable control law concept utilizing the Conmiged 
mixer as described above. This choice was made based on the 


following assumptions given in Reference 1 for implementation 
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feamehemnrecentrleturation control law process in an existing 
airframe: 
1) "The effector complement provides redundant effector 
systems and surplus control power for each flight 


Critical control force and moment." 


mee inert ltent control is a full authority, fly—by-wire 
ict pale tligit cantrol system.” 


Bee A control law design exists which has been carried 
out for the unimpaired airplane, but which is suf- 
Ptetentlyerebust that Only first order impairment- 
induced stability derivative changes need be 
Aecoumrted for in a drop-in reconfigurable mixer." 

Based upon these assumptions it is felt that the F/A-18 is 
well suited for the reconfigurablecontrol law study. 

In this thesis a linearized mathematical model was de- 
veloped which simulates the dynamic response of the F/A-18 
meceratt to stick and rudder inputs. The model includes the 
Memerol laws for the longitudinal, lateral, and directional 
axes for the cruise phase of flight. The fly-~by-wire actua- 
tors and sensor ee ee modeled and joined to the air- 
frame linearized, small perturbation model. The perturbation 
model was obtained from the Simulation Control and Technology 
meeoup, Plight Systems Branch, Strike Aircraft Test Direc- 
Mmenate, at the Naval Air Test Center. <A computer program 
was developed which, with existing software at the Naval 
meerestaduate School, will simulate F/A-18 longitudinal, lat- 
eel, and directional response to stick and rudder faputs. 

The model was designed to implement the control mixer 


Men Concept outlined above. Future work at NPS will utilize 


the simulation program to develop algorithms for determining 


13 


the control mixer gain matrix. In 
be used to study modern control augmentati 


Craft stability and geome oe 


14 


i ODE MeElHODOLOGY 


we ENT RODUCTION 

This chapter discusses the methods and assumptions used 
memtormulate a mathematical model which simulates the dynam- 
memresponse of the F/A-18 aircraft. The flight control sys- 
Me mweais described, including the simplifying assumptions used 
to develop the control system model. An overview of the com- 
plete system, which couples the control system model with 
mie airtrame small perturbation model, is then given with a 


Meet description of each functional component. 


Peer LiGHT CONTROD SYSTEM DESCRIPTION 

Peder aibededesccrapetomeor the P/A-18 control system and 
Mmmeory Of operation can be found in the Flight Control Sys- 
tem Design Report by McDonnell Aircraft Company [Ref. 3]. 
Mite following discussion briefly describes the basics of the 
Flight control system and the control law mechanization, and 
is intended to facilitate understanding the model development. 

The primary flight control system in the F/A-18 is a 
fly-by-wire, full authority, control augmentation system. 
The control law computations are performed by four flight 
@emerol computers operating in parallel. Each computer re- 
Memves input from the aircraft motion sensors, air data 
Mmomputer, and pilot stick commands. The computer operates 


feeeene input signals according to the control law algorithms 


IS 


and Bete the command signals to fly-by-wire electrohy- 
draulic servoactuators. Figure 2.1, taken from Reference 3, 
shows a functional block diagram of the flight control] )Syeieeue 
Exclusive of angle of attack and air data sensors, the system 
has quadraplex redundancy. The system provides two fail op- 
erate performance for augmented motion feedback control. A 
third failure causes the system to revert to either open loop 
direct electrical link control, or sStabilator meechanweem 
GOt tee Ole 

The control augmentation system is gain scheduled with 
angle of attack and air data to provide optimum flying quali- 
ties throughout the flight envelope. Cross axis intercon- 
nects (e.g., rolling surface to rudder interconnect) are 
provided for turn coordination and maneuverability at hag 
angles of attack. The control system also provides feedback 
Co eGoUunter anerrral coupling at high roll jrarece 

The F/A-18 has ten primary flight control Suptagecs 
Right and ett Stabilators, leading edge flaps, trailing 
edge flaps, ailerons, and rudders. Longitudinal contpoucee 
provided by collective stabilator, and collective leaiae 
and trailing edge flaps. Lateral-directional Compote, 
provided by differential stabilator, differential Meade 
and trailing edge flaps, ailerons, anderiddeqge Collective 
leading and trailing edge flap deflections are scheduled by 
the control laws and are a function of angle of attack. The 


flap positions are designed to provide optimum L/D during 


ale 
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S = Servo 
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Ps = Static pressure 

Qe =Dynamic pressure 


() = Trim motor 


—Pppe = indicates level of electrical redundancy 
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ae 
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True Airspeed, Ambient Temperature, 
Attitude, Heading, Selected Heading, 
Altitude, Altitude Rate, Vertical 
Acceleration. Data Link Commands 


Functional Block Diagram of Flight Control System 


cruise, and improve flying characteristics in maneuvering 


and high angle of attack flight. Figure 2°72) copied eee 


Reference 3, shows the control surface positions and direc- 


tion for positive deflect wom 


The objective of the thesis did not require a complete 


model of the F/A-18 control system as given in Reference 3. 


The following assumptions were made to reduce the model 


complexity. 


1) 


2) 


29) 


4) 


3) 


6) 
7) 


8 ) 


The aircraft is operating in the up and away flight 
phase. Under this assumption the control laws are 

Operating in the auto-flaps-up configuration. Con- 
trol law configurations for the takeoff amd) Pamcueee 
phases were not modeled. 


Only inner loop control is modeled. For inner oem 
control the pilot provides commands to the system 
The auto functions (outer loop control) were not 
modeled. 


Control is provided by the control augmentation sysimome 
The unaugmented modes such as direct electrical link 
or mechanical backup are not modeled. 


The failure logic provided to reconfigure the contimem 
laws in the event of a SenSOr Or ACtUatoOr |. a eee 
Noe Modeled. 


The aircraft is operating with gear up, speedbrakes 
in, and no -ext ernadgesiee eece 


Spin mode control logic is) not medeoled. 

The aircraft trim system iS not modewode 

High angle of attack conditions are not considered in 
the thesis model, therefore control law configurations 
for this flight condition are not modeled. Pore 


thesis high angle of attack is defined as flignhe@eome 
ditions above 15 degrees. 


18 


UOTZIITIIG BATATSOg Jo 

"spre AQTTIQUIS Ssarougap 

%( +) pure spxe Kpoq sazouop He) 
SUMONUS 910 SPUIPILTJJIOI Osae pur 
SUOPIIITJIP YILVFANS [oI WUOD AAPI;Soy 


ee 
CLC D RSH A Rea olP 


3p 


4sp 


UOTAOSATC pue suOoTy 


(1)  =SL1LON 


WY) 


a 


t 
. S 


LeOd OUTINGS TOTO] 


dis dior A 


Ze eigen 


4 


te 


iL. 


C. F/A-18 DYNAMIC MODED VG) Faia 

The aac tomal block diagram of the F/A-18 model which 
couples the flight control system to the basic airframe is 
shown in Figure 2.3. Thediagram represents a multi-input 
multi-output, sampled data, closed loop control system. 
Theory on the analysis of sampled data systems, as the one 
shown oa Fig. 2.3, is extensive and covered in 4 numberwos 
texts (see references). The development of the F/A-18 model 
assumes the reader has a rudimentary understanding of control 
theory, and in particular the theory of sampled data control 
systems. 

In the nomenclature used to represent the control signals 
pilot inputs are prefixed by the letter ‘P’. Actuating eae 
nals, and signals produced by the aircraft sensors are pre- 
fixed by the letter ‘'E'. Control surface deflections are 
prefixed by the letter ‘D'. The nomenclature used to repre- 
sent the aircraft perturbed motion variables and comerem 
surfaces is given in Table 2.1. To denote the motion axis 
which is being controlled, the signal will be Suftfixediiyagee 
either x, y, or z to denote longitudinal, lateral, oreceeeeee 
tional axis respectfully. Finally a matrix or a vector 
be denoted by an upper case letter. A scalar will be denoted 


by a lower case letter. 


The input vector shown ier ieee ee 


P(t)" = [px(t) py Cepia el (2 ai 
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PPS. 7/1 sO Ue rbetd yeoTg TeuCEscung ¢€72 Sanbig 


ca! 


(J) SXH=(3)Sq 


(A)A (4) Sq (4) Sq 


(9) 49+ (9) Sxd=(3)sx 
SOIWWNAG YOSNAS q/w 


(3)dd0+(3)WXH=(9)4 H 


(3)ea} O fxyeg (4) Og 
Z 


(3) ®XH=(4)d (4) Nd+ (4) OxXd=(H) Oa 
(A}A (3) a (4) d (4)d 
(3) d5+ (3) wx a= (2) wx (9) PAO+ (3) OXA=(9) eX 


(4) OG+ (4%) OxV=(T+H) OX 
SOIWVYNAG SWWeSUy SOIWYNAG YOLYNLOV 


v/a NI¥D SMW1 IOULNOD a/w 


Ar ee 


MOTION VARIABLES AND CONTROL SURFACERNOME TSS ieee 


q pi be hewate DEG/SEC 
nz normal acceleration G 

aa angle of attack DEG 
yr yaw rate DEG, SEC 
jee roll rate DEGys Be 
ny lateral acceleration G 
Sie right stabi Waimer. DEG 
Sic il left stabilator ‘J 
ter right trailing edge flap :: 
tel left trailing ved mea van : 


ler right leading edge flap 
el left leading edge flap 


Ue right aakere. : 
al ett 22 Wome a q 
ie 16 right rudder ig 
veal lett rudd er be 


represents the aircraft longitudinal and lateral stick, and 


rudder deflection in inches. The output ecerene 
E 
Y(t) “= [q(t) nz(t) aa(t) yr(t) rr(t) my(t) | eee 


represents the perturbed motion of the aircraft about some 
steady state operating condition. The motion variable Wage 
are degrees, degrees/sec, and G's. Each block in the dia- 
eram contains a mathematical model which simulates the dynam-— 
ics of that particular component. The controd a iaaae 
contains the aircraft flight control law algorithms modeweg 
as linear, time invariant, discrete state €quae moms. One 
processing channel of the flight control computer describes 
above is represented in the control law model. The input 
vectors to the control law model are the discrete Stickam 


rudder input signals in inches, 
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c \ 
Pee = ipxCk)  py(k) pz(k)] ea) 
and the discrete motion feedback signals from the sensors, 


ately = | a Cemua@emaacGk) yr(kjerr(k) ny(k) ] Cre | 


(The units of the motion feedback signals are degrees, degrees/ 
sec, and G's). The input vector to the control law equations 


is therefore 
e 
Go) ee REeGr) mor Ce) e) 
Mime Output vector from the control law block, 


Ee(k)° = [estr(k) estl(k) eler(€k) elel(k) eter(k) 


etel(k) ear(k) eal(k) err(k) rama) 12.6) 


mepresents the discrete command signals to the flight con- 
trol actuators in degrees. The actuator penne nd Signals en- 
ter the GAIN block which represents the configuration gain 
matrix discussed in the introduction. The state space equa- 


mmons Shown in the actuator block model the dynamics of the 


meme t control actuators. The input vector to the actuator 
Block, 
IE 
Pum esctror) esti(t) e@ler(t) elel(t) eter(t) 
SeeieGtimear( th eal(t) err(t) erltt ) } ee. 7.) 


Mepresents the continuous time, actuator command signals in 


Meemees. The output vector from the actuator block 
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represents the control surface detlections tie@eoreee 


D(t)” = [dstr(t) dstl(t) dler(t) dlel(t) dter(t) 


dtel(t) dar(t) dal(t) dru ean ae (27am 


The deflection vector is input to the airframe small pertur- 
bation model represented by the state variable equations in 
the airframe block. The output vector from the small per- 
turbation model, Y(t), enters the sensor dynamics block ieee 
contains the state variable model for the aircraft rate gyros, 
accelerometers, and angle of attack sensors. As discussed 
above the sensors output the feedback signals which are sent 
to the control laws. In Chapter III the mathematical models 
which simulate the dynamics of each component in Fig. 203 

will be developed in detail. 


Analog to digital converters are modeled as impulse sam= 


poker s+ It is assumed that all samplers are operating at the 
Same, constant sampling rate. (The actual system uses multi- 
rate sampling. In the thesis model only a Single Samples 
rate is used. The program actually allows any desired sam- 


pling rate to be input.) The mathematical operation per- 
formed by the impulse samplers is shown in Figure 2.4. 
Modeling the analog to digital converters as impulse samplers 
is a valid assumption if the quantization error of the ace 
system is at an acceptable level. Digital to analog conver- 
ters are modeled as zero order hold devices. The mathemati- 


cal operation performed by the zero order hold is shown in 
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Maes oe fF inaliy if is asstmed that the processing time be- 
tween the sampler inputs and zero order hold output is very 
much less than the sampling period (i.e., no processing delay 
time is assumed). This assumption is used in transforming 
the continuous time state equations (actuators, airframe, and 


sensors) into discrete time equations. 
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f 


Figure 2.4 Mathematical Operation Performed by Impulse Sampler 


eee ee 
oo ee eee ee 
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Figure 2.5 Mathematical Operation Performed by Zero Order Home 
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ere oO UCbeo Ey ELOP MENT 


Percoll ROL G@AW MODEL) DEVELOPMENT 

imine Midtaaiebareral, andydirectional control law 
models were developed from the block diagrams and information 
provided in Reference 3. In addition to the simplifying as- 
eimptions listed in Chapter II, many of the ea onehit ial 
this system were eliminated based on the following 


memsiderations: 


1) The control law model will be coupled to the small 
perturbation model of the F/A-18. It is assumed 
Enameoperataonm of the control system will remain 
within a linear region. Therefore the non-linear 
components in the control laws were eliminated. 
treseminchiwde position limiters, rate limiters, 
and dead band regions. Other non-linear functions 
in the system which were essential to the model 
(pitch stick gradient for example) were linearized 
by a taylor series expansion. 


2) For the same reasons discussed above, the portions 
of the system which provide inertial decoupling 
were not modeled. 


3) In the model, the control signals are input as dis- 
crete signals. Therefore stick and rudder dynamics 
are not modeled. 


4) Noise, which may be detrimental to control system 
Eetitormance due to aliasing, iS not introduced into 
the model. Therefore the anti-aliasing prefilters 
were not included in the model. 


pees ihe structural modes were not included in the F/A-18 
rebar rame mode | . iMfjeuwecihonre. Ehewstructural notch») ti1-— 
BetiesSeWere Tot included in the model. 


fee fO prevent discontinuities in the Signals, the con- 
Seow avomleblize taders in portions of the system. 
PiscomemmuinteLes Gould O¢cur during start up, fail- 
ecm or Eranusitions. since none of these condi- 
tions are included in the thesis model, faders have 
not been modeled. | 
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7) The F/A-18 control system uses multi-rate sampling 
in the input and feedback paths (20,540) andmoge 
sampling rates are used). To develop a state var- 
liable model only a single sampling rate was con- 
sidered. Therefore the iteration averagers, used 
to mathematically combine two discrete signals of 
different sampling periods were not modeled. The 
Simulation program allows any desired sampling rate 
to be input. For this thesis 80 hz was Used eae. 
Sampling seat ee 
Figures 3.1 and 3.2 show the simplified block diagrame 
of the longitudinal and lateral-directional contro lm 
The inputs are the discrete stick and rudder commands, and 
motion feedback signals. The output signals are the commands 
to be sent to the flight control actuators via the reconfig= 
urable gain matrix. Together Figs. 3.1 and 3.2 make up the 
control’ lawtblock “show mine ton. ee 
The blocks in the control law model represent two basic 
transfer functions: Function gains, and digital tata 
Table 3.1 lists the notation used to represent the transfer 


functions. Function gains and digital filters are descmieg 


below. 


TAB EES Set 


CONTROL LAW TRANSFER FUNCTION=NOTs ies 


Page i lex Transfer ftmetaon 
oe Punetion «cam 
Eee Longitudinal digital? ee 
R Lateral digital ei igitees 
a Directional digitale: wee 


ass 


Hong? eudinal 


biek coe 
2, Sere gra ets 
Ss 
Longitudinal Gomlective 
tabilator 


Input 


Normal 


meceleration 
pe 
fe 


q(k) , 
‘ ie 
| po F68 | F32a j 
Patch Rate qc qc or 


———= 


Collective 
aa(k) Leading Edge 
P rap 


Angle of 


Attack CeGwee nap Cicacd 


etex(k) 


~405 (a a 


Collective 
ira rngvedge 
Elap 


Braves emo limpiered Longitudinal Control Law Diagram 
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plateray. 
EV LCE Jas \ aloe 
py (k) aE, SZ ae . FY 3 eleva 
ae el ae : 


Lateral Input] ~~ ri,ps,nz| Leading 


Diftfereneiaw 
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eteya 
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Bi esty (k) 
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rk 
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rrik) 


roo 
Pp ic, Ps, 
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Aelia ae 
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Figure 3.2 Simplified Lateral-Directional Control Law Diagram 
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Pare une ELOns 

Mipewruneti ons imetine control law diagrams perform the 
feietem gain scheduling described in Sec. II.B. The functions 
operate on the air data, angle of attack, and normal acceler- 
ation to compute the system gains. Table 3.2 lists the nota- 
tion used to represent the function inputs. The mathematical 
equations which define ore Poepooredubesvare given an Refer- 
ence 3. Copies of the functions used in the control law 
model are given in Reeney aioe srmularron program the 
Mmimetilon gains are computed using steady state conditions for 


all input values. 


Ni 8 1b MS Re 


PUM ELON INPUT VARIABLES 


ps Pidtieated Static pressure ey ob t2 
qc Dynamic pressure Meise stot 2 
ri Pressure ratio (ps /qc) ND 

nz Normal acceleration G 
alpha Angle of attack DEG 


2. Dieiiea Pekilters 


Lead-lag filters are used in the system to shape the 
output response and provide nacre gain and phase margins. 
Meintegrator is used in the forward loop of the longitudinal- 
System to provide zero steady state error between command and 
Meeadback. The control system design report [{Ref. 3] gives 
Mmeetilter’s continuous time transfer function. For the con- 
Meer law model, the digital filter coefficients were computed 


using the Tustin transform. All filters were modeled as first 


od 


order systems with two numerator coefficients and a single 


denominator coefficient. For example 
H(Z) = PONT 7 (37 
Z = Pop 
is the model for filter P9. Appendix B gives the dioieam 


filters used in- the model, and the method used to-compute the 
Ccoetticientrs. 
3: otate Space Modells 
Standard control system analysis techniques were used 
to derive the state space models for longitudinal and lateral- 
directional control laws shown in Figs. 3.1 .amg3 2 ee 
following sections outline the procedures used. 
a. Longitudinal System 
The block diagram in Fig, 3.1) contain eas 
input/output signal flow paths: 
1) ~-Pitch rate to collective stabwla.. 
2) Normal acceleration to collect: vemem-ipeel-neor 
3) Longitudinal stick to collective stabi, 
4) Angle of attack to collective leading edge flap 
5) “Ang le on eae backer te Ne vieticter trailing edge flap 
To obtain the individual path transfer fumet ions meme signee 
are mathematically combined to give the following three Z- 


Lranslorm edvuettogs. 


Estx(Z) = H1(Z)Eq(Z) — H2(Z)Enz(Z) - H3(Z) Epx (Ze 


Elex(Z) = H4( 2be ee (3a 


peers) = lst Z)haaCZ ) (3.4) 


The expressions for each transfer function are given in Ap- 
pendix C. To obtain a state space expression for longitudi- 
meercontrol laws, the individual Z-transfer functions are 
first expressed in state variable form. The state variable 
equations are then combined according to equations 3.2 - 3.4. 
Mis procedure is as follows: 


Mme livewttemerter function for the pitch rate to collective 
Stabilator path can be expressed as: 


(eee WOES uae Cosy 


(Z-P9D) (Z-P2D) 
The numerator coefficients are functions of the system gains 
and filter coefficients. The roots in tne denominator are 
the poles from Piiteemcoug@and) hf 2. «Appendix C eaves the de- 
tailed expressions for the numerator coefficients). 
2) The state space representation of Eq. 3.5 is obtained 


using the parallel programming method outlined in 
Reference 4. 


xl(k+1) ON 10 ke) il 

-| il * Jexce (eG) 
eke OMe 2 Wie 2 Ok) 1 
ell) 


nes) fqsed ese al 


|: qst3 eq(k) Ost) 
ZK } 


Applying the same procedures to the remaining Z-transfer 


Memeetons in Bq. 3.2 results in similar expressions: 


rs) 


x3(k+1) POD O Aika i 

|: | J-[ Jerse (3.8) 
x4(k42) =. 0 PPS a 1 

xieG@i)) 


estx2Z(k) anes el | |: nzst3 enzce (328 


K 2 le) 
for the normal acceleration paths i204) 


x5(k+1) = POD x5(k) + 1 px(k) (3.10) 


estxSGk) = pxstl x5(k) + pxseZ paw (3 = 


for the longitudinal Stick path ice 
Appendix C details the procedures used to com- 


pute the coefficients in the output equations. The nomena 


ture used to represent the coefficients in the output equation 


combines the notation of the input Signal an@eemepme 
control surface, followed by a number indicating the coeffi- 


cient's numerical order in the equation. For example: 
GStl, @s &2.00G cao 


are the first, second, and third coefficients in the page 
rate to stabilator output equation (Eq. 3.57 Wate 
system of nomenclature the respective signal path of the co- 
efficient is easily identified. 

3) The state variable equations for H1(Z), H2(Z) andi 


are now combined according )to) Eanes 


estx(k) = estxl(k) = estx2(k ))— ecmaree (3.0 
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meee Eqs. 3.7, 3.9 and 3.11 gives 


mek) 
x2(k) 
estx(k) =[qstl qsi2onzstl=mzst 2 epesie | |ecc (kc) + 
coor kk } 
Con ep 
eq(k) 
fased-nzeeslf |+ -pse1 px(k) Ger 3) 
| enz(k) 


mer the output equation. RQppeeonms oO, o.0,.-and 3.10 can 


be combined to give 


x1l(k+1) Pop pe Ck) Lae O 

x2(k+1) Ee sok ) 1 Oteq(k) O 

mo(k+1){= POD x3(k) +10 +]o]px(k) 
x4(k+1) Pat we ke) Oelilenz¢(k) fo 

x5(k+1) PODHW x5(k) O 0 O Cele) 


The state space representation of equation 3.2 is given by 
Mes. 3.13 and 3.14. These equations output the collective 
Peeabilator command given the longitudinal stick and motion 
feedback inputs. Note that the motion feedback and longitu- 
Menal inputs have been separated. This facilitates coupling 
imme control law equations to the aircraft equations to be 
developed later. 

Mmepeoitilar State variable equations are derived for the 
angle of attack to collective flap path transfer 
mincomons Haq 2). and H5(Z): 

Morir =P lID xo(k jo + 1 -eaa(k) (on. 5:) 


emer ljw=saalel=x6(k) + aale2 eaa(k) Goo.) 
for the AOA to collective leading edge path H4(Z); and 
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| 


x7(k+1) = P12D x7(Cejeee te eeenee (3. 1s 


etex(k) = aatel x/(k) + aate2 eaark) (3. ee 


for the AOA to collective trailing edge path H5(Z). 


5) The state equations for collective stabilator, sama 
collective leading and trailing edge flaps are now 
combined to give the following state variable model 
for the longitudinalwycone ve ae. =. 


x1l(k+1) POD ral) 
S2Cke is P2D x2(k) 
x3(k+1) POD x3(k) 
x4(Ck+ bia P5D . x4(k)] + 
Yc POD x5(k) 
xO(k+1) Pavia HWx6(k) 
x7(k+1)] P12Di}x7(k)| 


im Ome C 
EO. 0 
O 1 Oleq(k) 0 
O 1 Ojfenz(k) +]O0]px(k) (3.19) 
O 0 Oj] eaa(k) ip 
Dp. © 1 0 
© Of 0 | 
x LG 
x 2 ke 
estx(k) gstl “qst2 —nzstil=n2zst2 —pxs eee 0 x3 Cra 
viex(i)]-{ O 0 O 0 0 aalel 0 eaes - 
etex(k) 0 0 O 0 0 O aatell[/x5(k) 
x6(k) 
x7 Gis Wy 
qst3-nzst3 O eq(k) —-pxst2 
0 O etl ou ene (i) 0 px ke (3.205 
O 0 aate2ileaa(k) O 


Equations 3.19 and 3.20 are the state space representation 
of the longitudinal control laws shown in the block diac 
in Fig. 3.1. The equations represent the longitudina leome 


trol law model which computes the collective stabmeaaes 
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Sonmand, and collective flap commands, given the longitudinal 
Stick and motion feedback inputs. 
b. Lateral-Directional System 

The procedures outlined above are applied to the 
mm 3.2 to obtain the state space model for the lateral- 
directional system. (To complete the discussion of the con- 
trol law model these procedures will be briefly described. ) 
iite individual path transfer functions are first a uies med by 
fienematically combining the signals in Fig. 3.2 to give the 


mMotlowing equations: 


Esty(Z) = -H6(Z)Err(Z) + H7(Z)Py(Z) + H8(Z)Pz2(Z) (3.21) 
Eley(Z) = -H9(Z)Err(Z) + H10(Z)Py(z) Gree) 
Etey(Z) = -H11(Z)Err(Z) + H12(Z)Py(z) (cae) 

Ea(z) = -H13(Z)Err(Z) + H14(Z)Py(Z) + H15(Z)Pz(Z) (3.24) 


cio Ay ey a ( 2) are Hl/(Z)Err(Z) + 


PUet@Ze mae? ee O( 2) Py (4201 Z)Pz(Z) (3.25) 


These transfer functions represent the following input/output 
Signal. paths, numbered respectfully: 

Meeeeroll rate to differential stabilator 

Mmieateratestick to ditterential stabilator 

pee Rudder pedal to differential stabilator 

©) Roll rate to differential leading edge flap 

my Lateral stick to differential leading edge flap 

Meee koll rate to differential trailing edge flap 


Mee Lateral stick to differential: trailing edge flap 


ey 


13) Roll rate’ to arlewon 

14) Lateral stick to aileron 

15) Rudder pedal to aitewen 

lo) Yawerates te ridden 

1/7) Roll Yate tor ridamer 

18) Lateral acceleration to rudder 
19) Lateral stick@tovreudde: 


20) Rudder pedal to rudder 


The expression for each transfer function and corresponding 
State equation are given in Appendix C. Note that the dif- 
ferential stabilator, ailerons, and rudder contain the Ga 
fer functions for the cross axis interconnects (é€.¢9., Honea 
and H15(Z) represent the rudder to rolling surface intemeemes 
nect, and H19(Z) represents the rolling surface to rudder 
INberconneck..) 

The state equations for the lateral and direction- 
al control laws are given as: 
esty(k) O -rrst 
eleyCk)\)s=10 -=rrite 


etey(k) Oi 2 rite 
ca(k) QO -rra 


pyst pzst 
eyr(k) +}pyle 0 pvy(k) (3.268 
err(k) pyte O pz(k) 
eny(k) pya pza 


ae) 


for the lateral systemiyamnd 


x8(k+1) Y3D K8(k) 7 
x9(k+1) NOSND, x File) 

xlOCk+1) Y5D SONG c 
X11 Ck Yo) po lige 
xl2(k+1) Y5DNW& 2 Ck 
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iene c Sey 0 
Cuero Neyr Cte CO s0aipy Ck) 
O) stl 10h | Bseie Ole) 0 0 O3e 27.) 
O O Offeny(k)J+]1 O]pzck) 
0 O O O 1 
co) 
x9(k) 
er(k) ={-yrrl aoe leer rr 2 oy rl pzr-] slbOik ) 1) + 
ol Oke) 
Sek) 
ENG ol, py(k) 
(Sie ages ey dfersce) | [pyr2 pzrz] Co 28)) 
eny(k) pz(k) © 
for the directional system. Note that the lateral system is 
Bemorder zero. No filters were included in the model. 
wn eombanine the Control Law Models 


hetimomso 7 aud 3.270, and 3.26 — 3.28 are com- 


bined to give the 3-axis control law model. 


eq(k) 
enz(k) 
Zsa Mex 2 2x1 12x6 eaa(k) oe) 
Metelc 1). = heme Celia) se bwemeyaCk) = be py k) (3.29) 

erm (ic) pz(k) 
eny(k) 

estx(k) eq(k) 

elex(k) enz(k) 

etex(k) Sele le x 1! exo eqatCk ) Sop x (ie) 

Sere «k) = Ce Xce(k) + Dfe eyr(k) + Be py(k) Cores 0) 

eley(k) err(k) pz(k) 

etey(k) eny(k) 

ea(k) 

er(k) 


mentions 3.29 and 3.30 are written in terms of the matrix 


meer icients in Appendix D. These equations are represented 


a 


by the discrete state equations in the control law bloc aa 
PUINCe wee on: 
3. Variable Gain Manrix 
As a final step in the control law model development, 
the variable gain matrix is introduced. The command Signals 
which are output from the control laws (Eqns. 3.29 and) see 


are distributed to the right and left actuators accordimeuee 


the diagram in Fig. 3.3. For the unimpaired aircerat tee 
individual gains in Fig. 3.3 are set to unity. - The fou 
ing matrix gain equation represents the diagram in Fig. 3.3. 


'GAIN' 

eskric. 1 0 0 =I 0 0 O 0 

estl(k) 1 0 0 1 O 0 0 Otlestx(k) 
eler(k) 0 il O 0 il 0 0 Ollelex(k) 
elel(k) O 1 O Oe O O Oleteena 
Sten kt —1ee 0 1 0 QO -1 0 Oflesty(k) 
ete File) O O at O O ll O Olleley(k) 
ear(k) 0 O 0 0 O QO -1l O}letey(k) 
eal(k) 0 O O O O O It O llea(k) 
err(k) 0 O 0 O O O O 1 ler (Ck) 
erick) O O O O O O O Jk 


(3 303 
The matrix gain equation will be recomputed for the impaired 


aa lead ee 


B. ALKCRAFT MODELODE VELLORE a: 

A 3-axis control law model has been developed which will 
Operate on the stick and rudder inputs and the motion feed- 
back signals. The model outputs the control signals which 
are distributed to the aircraft actuators through the vane 


ble gain matrix. The next step in building the F/A-18 model 
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Figure 3.3 Control Law Command Signal Distribution and Gains 
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is to develop the state ja Puan equations for the air? rane 
actuators, and sensors. 
1. Airframe Model 

The airframe model, obtained from the Flight Systems 
Branch at NATC, is a linearized small Pon model for 
both the longitudinal and lateral-directional modes. The 
dynamic stability and control derivatives were generated from 
NATO 6 F/ Aas Sanna ies Package. 


The aircraft state fGiaeeee 


model is trimmed, unaccelerated l-g flight. The airframe 


State equations are given as: 


LONGITUDINAL 
u(t) 4x4 u(t) 4x3 dstx(t) 
w(t) = Fx w(t) + Gx. dilexis (3.3 
q(t) “Cie, dtex(t) 
thed(t) thed(t) 
aes, Cee Gee) 3x3 dstx(t) 
noe) = Hx w(t) + Ux dle Ga (3.338) 
aac q(t) dtex Cen 
thed(t) 
LATERAL DIREC! LOM 

dsty(t) 
aE) eC oe) 4x5 dley(t) 
gn Oe ae ge gy A ie) + Gyz dtey(t) (3.349 
pte) pe» dar 
ellgtal1 Ge), Dice alien. (ie) 

dst 7 Gee 
Ee) Bee ae 3x5 dley(t) 
p(t) = Jie Ze Gre + Dyz “deevan (3.338) 
ny(t) Bae) da(t) 
Pitter stem) dr(t) 


Airframe variable definitions and units are Ji steq. 7 eee 


3.3. The definitions of the stability and control deriva 
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(ei 323 


PRP RATE VARIABLE DEFINITIONS 


hotation Variable Uinieine Ss 

w(t ) longitudinal velocity perturbation eyes 
w(t) perturbed normal velocity fey 6S 
q(t) perturbed pitch rate Ead7/s 
Ce ) perturbed lateral velocity ft/s 
yr(t) aircraft perturbed yaw rate Pay) S 
ao t) aircraft perturbed roll rate rad/s 
thed(t) Jiserartepereurbed pitch angle rad 
phi(t) Meeerart perturbed roll atele ‘ 
etx ( t ) Sclmeceuvo sStabiuekatror deflection i 
dlex(t) eollective leading edge flap deflection " 
dtex(t) collective trailing edge flap deflection " 
dsty(t) @ifterential stabtlator deflection - 
dley(t) (iirerential le@adums edoe deflection - 
dtey(t) differential trailing edge deflection ‘“ 
da(t) elem omeder Vecrsa on : 

fier ( t ) elilalepeie (ao edevenemliong i 

ma ( t ) Memmalemaece eration fay SZ 
= 6 (oe angle of attack rad 
ny(t) lateraleacce eration egay, 7 


and associated units, which make up the matricies in Eqs. 
Sees through 3.35 are given in Appendix E. Note that the 
units of the perturbation model are not compatible with the 
Meemators or sensors. (The actuator output units are in 
degrees, and the sensor input units are in degrees/sec, de- 
ees, and G's.) The input and output variables for the 
Perturbation model were scaled in the simulation program to 
Meeperly interface the models. Combining equations 3.32 - 


Peo. gives 
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poe ico) u(t) dstx 
w(t) 4x4 4x4 w(t) 4x3 ax 5) ed ex 
att) F x O Greats) Gx @) dtex 
thed(t) = --------- thed(t) + --------- dsty (3. Ste 
v(t) Ax4 4x4 v(t) 4x3 4x5 dley 
ig G ie 9 0 BY Zones 0 Gyz dtey 
p(t) pGe da 
epotat Gee, eum Ge), dr 
(Gee doa 
Crs 3x4 3x4 wr SS 3x5 dlex 
mZee.) Hx O a Ie Dx O diex 
aa(t) = --------- thed(t) + --------- dsty (3.30 
Et i bee 334 Aan) Cpe 3x5. 4d Dew 
DCe O hy 2. Vay O Dyz dtey 
my eo ie ee) da 
oval (Cc, dr 


To study reconfigurable flight controls the aipenmaae 
model should be capable of using the full set of control sur- 
faces available to produce the required forces and moments. 
To achieve this the control surfaces are split into indepeme 
dent right and left hand complements (i.e., right elevator, 
left elevator, etc.). The equations are then coupled so that 
a complement of control surfaces used either collectively, 
differentially, or as a single side, will produce the appwge 
priate moments. For example the stabilators deflected col- 
lectively will produce a pitching moment, deflected 
differentially will produce primarily a rolling moment, “anime 
Single side deflected will produce, to some degree, moments 
about all three axes. 

The F/A-18 airframe modeled in equations 3.30 9am 


3.37 inherently offers control coupling through the Stabe 


tor, leading edge flap, and trailing edge tlap suriaces™ 
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Mi icional Pit aime omp line could “bewachieved with the 
meberons and rudders. At one Cin eee was written 
Monmtrol derivatives were not available on the longitudinal 
effects of the ailerons or rudders. (The rudder is capable 
ee toe-in or es one and will effect the longitudinal re- 
meemse Of the aircraft. this feature is normally used during 
takeoff and landing.) 

Poe phat omeotmnal Suttace derlections into right 
and left hand complements, the following equations are used 
[meen Compute the deflection inputs to Eqs. 3.36 and 3.37. 
Meso refer to Fig. 2.2 which shows the control surface posi- 


tions and corresponding positive deflections.) 


LONGITUDINAL DEPLECLLONS Poth weer lONAlD DEP LECTIONS 
Meex =. (dstl + dstr )/2 ee dis be ag Ss tT 
dlex = (dlel + dler )/2 dley = -dlel + dler 
mmex = (dtel + dter )/2 Geeovw —sateleo= -dter 

ree edie merely ) / 2 

pee da eee ) / 2 
Meere t and 1 correspond to right and left surfaces. These 


Meret Ons are rewritten in the following matrix format: 
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| 


dstx 
dlex 
dtex 


dsty 
dley 
dtey 
da 
ia 


| 


Ore Ca 


| 


ES GS 0 HS 


© Qo 


Oo oor 


ae ES ae 


t 


d1S tebe ron 


EO Nee Miao tae 
0 O O 0 O O 

5 O O O O O 

Ci a5 5 O O O 
LATD ‘distribu ticumnap aie 
0 O 0 O O O 

=i O O O 0 O 

QO =i 1 O O O 

O O QO: Sh i O 

O O O O ee 


eG Cnr) 


G7 ES Se TE Se FE. 


die teas 
dst 
dler 
dlel 
dter 
dtel 
dar 
dal 
alle gar 


alg ewil 


dstr 
dsl 
dler 
dlel 
dter 
dtel 
dar 
dal 
dice 
det 


(3.389 


(3.339 


Note these equations are for the unimpaired aircraft only! 


Damage to one or more of the control surfaces will change the 


LONG and LATD matricies directly. 


Replacing 


Wi thaterme so. hese 


of equations 3.34 .amdeee 


the input vectors in 


ing modified airframe equations: 


thed (n>) 


yD) 
Pee.) 
p(t) 


Dilek) 


Fim 
4x4 | 4x4 
Fee, 370 
4x4 | 4x4 
Oe yz 


u(t) 

w(t) 4x3 
q(t) Gx 
thed(t) + 
v(t) 4x5 
yr(t) Gyz 
Tae nie) 

DhaeCe) 
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oh 


EGse 


3.36 and Ga 


39 gives the follow- 


dstr 
ds teal! 
dler 
dlel 
dter 
dtel 
dar 
dal 
‘alg arse 
drill 


(3.405 


Hm Dm alshese 


u(t) al ied 
at ) 3x4 | 3x4 w(t) oo Od hem 
mz t ) Hx | O ea, Dx LONG dlel 
aa(t) = --------- thed(t) + --------- dter 
met) 3x4 | 3x4 v(t) 3x5 5x10 dtel 3.41) 
leet) O | Hy aye & } Dyz LATD dar 
my (t ) i a a deal) 
pines ) alieje 
aljaull 


meee names of the individual matricies appear above the 
equations. 
Pe Actuator Model 

Micon Uiekwonceamer the L[lisht control actuators 
were given in Reference 3 and are listed in Appendix F. 
These transfer functions are low order approximations of the 
more complicated actuator models presented in Reference 3. 
The models were used in the F/A-18 rigid body stability anal- 
ysis described in Reference 3, and approximate the frequency 
mmeracteristics of the higher order models out to 5 hz. 

[omic ouneareare Ene actuators imto the modified air-=- 
imeame model, the transfer functions are first put into state 
variable form. This procedure is outlined in Appendix F. 
The following equations represent the actuator state variable 


model: 


24x1 24x24 24x1 24x10 10x1l 


V@meie=o Fa kat). + “Ca Fat) (Come 2’) 
OB eal 2 & Dex |. 
Dit) = Ha joi Ginn) se.) 


si 


Ea(t)* = festr(t) estl(t) eler(t) elel(t) eter(t) etel(t) 


ear(t) eal(€t) errCe eee 


p(t) = [dstr(t) dstl(t) dler(t) dlel(t) dter(t) dtel(t) 


dar ete) dr Ee) driCe jal 


The input vector, Ea(t), represents the input signal in de- 
grees from the control laws, via the GAIN matrix. The outpue 
vector, D(t), represents the right and left control Suri 
deflections in degrees. Equations 3.40 - 3.43 are combined 


to give the following equations of the airframe plus actuator 


model: 
De NOC, Gp 
8x1 8x8 | 8x10 10x24 8x1 8x10 
Ct oe Fm | Gm Ha Xm(t) O 1O0xa 
—----- =  ----------------- ----- + ---- Ea(t) 
atl 24580 | ote Die 24x10 (3ooee 
Ceep, 0 | lp ee) Ga 
Hp Ooi ie ) 
q(t) 
OZ) 8xl 
aatr 6x8 | 6x12 12x24 Xm(t) 
ie): Hm | Dm Ha 9 ----- (3.4m 
RigGe”) 24 
Reve ) Xa(t) 
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Fs before the names of the matricies appear above the equa- 
rons . The arriranemplus actuator model inputs the command 
Menals £O tame actuators, and outputs aircraft motion. Since 
the control surface deflections in the perturbation model are 
in radians, it was necessary to scale the Hm and Dm matricies 
to interface with the actuator model which outputs deflections 
in degrees. This is done in the simulation program. 
fee oensor Model 

iigeene 1 ina! ees of the development of the aircraft 
model, the state variable model for the aircraft sensors is 
incorporated into the airframe plus actuator model. The 
sensor transfer functions are given in Appendix G. In state 


variable form the sensor equations are given as: 


gig) 
Ma) 
Lex xe eel ial i Ome at) C3446) 
XsCt) = F's Xs(t) + Gs et) 
p(t) 
ny(t) 
Sait ) 
enz (te 
eaa(t) le lela i 
eyr(t) = Hs KeCt) (32473 
epee) 
ete) 


Whe sensor model inputs the aircraft motion variables in 
degrees/second, degrees, or G's, and outputs the correspond- 
MmoeestisnalsS to the control laws in the same units. Combin- 
ing Eqs. 3.44 - 3.47 gives the following model of airframe 


plus actuators plus sensors: 
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Fps Xps(t) Gps 


Sp al 32x32 | 32x eee 32x10 
Xp(t) FP 0 Xp(t) Gp geal 
----- = -----------------  ----- + ----- Ea(t) (3.48) 
eee 11x6 6x32 | Si ele 11x10 
Xs(t) Gs Hp | Fp Xs(t) 0 
Hps Xps(t) 
eq(t) Sa 
enz(t) 6x32 | 6x1l Xp(t) 
eaacty = 0 | Hs --~-- (3.495 
ey PCr) 11xl 
Ore ce) 7 Nes (ie) 
eny(t) 


The aircraft model inputs the actuator signals from the weg 
trol laws via the GAIN matrix and outputs the motion signals 


from the sensors which are sent to the control law €quatasaees 


GC. ASSEMBLING THE VO ERA Lia sys EM CREE 

A mathematical model for each component in the control 
system shown in Fig. 2.3 has now been developed. Before the 
individual components of the model can be assembled, the dis- 
crete state equations for aircraft model (Eqs. 3.48 and 3.49) 
must be computed. Performing this operation the discrete 


State equations for the aircraft are given as: 


Ae Salk 4 364 3) Aa 43 x0 ele ae 
Xps(k+l) = Aps Xps(k) + Bps EaCk) (3.3508 
6x1 6x44 ap all 
Es(k) = Hps Xps(k) (3.4 
Where 
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Xps(k)* = [Xm(k)|Xa(k)]Xs(k)] 


Ea(k)* = [estr(k) estl(k) eler(k) elel(k) eter(k) etel(k) 


ear(k) eal(k) err(k) erl(k)] 


Meck)" = [eq(M) enz(k) eaa(h) eyr(k) err(k) eny(k)] 


The Aps and Bps discrete matricies are computed as follows: 


fps = a FPS es Coe 
ts st 
bps = f a BS “ds X Gps Corea S)) 


Where ts represents the svstem sampling time. 

The GAIN matrix is now introduced to interface the dis- 
crete aircraft equations with the control law equations. Re- 
Meeacineg the input vector, Ea(k), in the aircraft equation 


Meee the r.h.s. of the GAIN equation (Eq. 3.31) gives the 


moelowing: 
ox 1 Jn a ex Il ep LO 6c x | 
Xps(k+l) = Aps Xps(k) + pps  GalN Beck) Cl) 
ox | 6x44 dee | 
eck) = Hps Xps(k) (Gciatas 


For convenience the control law equations are repeated 


below: 
al le? ale b2xor Ox) ee 5 5 xe: 
meck+1) = Ne eG): 3 ieee Bic Es(k) + emer (eo BDO 
ex | cL eel Sx6. 6x1 8x3 Spell 
lewk) = Ce Xce(k) + Die test 1k) Dew eee 23.30) 


Sal 


Where: 


eth = “eccen eaae 

Es(k)° = [eq(h) enz(k) eaa(h) eyr(k) erecta 

P(k)~ = [px(k) py(k) p2(k)] 

Ec(k)" = festx(k) elex(k) etex(k) esty(k) eley(k) etey(k) 


eat keer om 


Equations 3.52, 3.53, 3.29, and 3.30 are OW com> emecmnes 


give the following matrix equation: 


43x1 43x43 43x10 10x8 8x6 6x43 ; 43x10 10x88 8x12 ae 
Xps(k+1) Aps + Bps GAIN Dfc Hps 1 Bps GAIN Cc Xps(k) 
nessa 12x6 6x43 | 122 1 2a 
Mec es ) Brew ios | Ac XcCe 
43x iO nO <6? oe 
Bos: —GADN) Dewr a 
$+ Re ---------- py(k) (3.3548 
12x38 pac) 
Bc 
6x1 6x43 | = 6x12 43x1 6x3 
Es (k) Hps | 0 Kos Ck) OC pei) 
~----- = ----------------- -~----- + --- py(k) (3.35@ 
8x1 8x6 6x43 | 8x12 12x1 8x3 pz(k) 
Beds) Dice ps | GG cok Dec 


These equations model the dynamic response of the F/A-18 


System shown an Figure 2 


D. MODELING EFFECTOR IMPe tien 
Effector impairment is divided into four groups tenes 
‘effector impairment classes’ (EIC) [Reft. 2) fb wee 3 a 


copied from Reference 1, defines the ELC S¥andung peer 


a2 


CONTROL STABILITY & 
ALEECIED SOMmtTK CL ArPFECTED 


ACTUATION LOSS, 
SURFACE FLOATING 


ACTUATION LOSS, 
SUREACc GEN ERED 


Ele =F 


PARTIAL SURFACE 
EOSs 


Pigure 3.4 Effector Impairment Classes 


a 


their affect on aircraft stability and COnETO ieee eee 
Fig. 3.4, effector impairment class one is unique in thee 
only aircraft control is affected. The remagnane in pene 
classes affect aircraft stability as well as Gone. come 
self repairing control system must be capable of detecting 
and classifying effector damage to compute the proper recon- 
Lilguracion Gains. 

For the case of EIC=1, one or more of the aircraft con- 
trol derivatives will be altered. In the F/A-18 model de- 
veloped above this is reflected in the LONG and LATD matricies 


(Eqs. 3.38 and 3.39). For example if the right stabi latoueeeee 


impaired the elements LONG(1,1) and LATD(1,1) would be seumige 


Zeros 


LONG 
QO 5 O O 0 QO Q 0 O QO r 
0 Os<5 5) 0 OQ 0 QO Q 0 
0 O OQ O Bee» 2) 0 0 QO 0 
LATD 
QO 1 0 O O 0 OQ 0 Q 0 
OQ O 1 -l 0 0 O 0 QO 0 
O 0 0 Ce el il OQ 0 0 Q 
O O O 0 O O -l 1 OQ O 
OQ QO 0 0 OQ 0 0 0 4. We 


Compare the above matricies with Eqs. 3.38 and 3.39 foreeme 
undamaged aircraft. Note that the damaged system will now 


produce a lateral input for a given longitudinal command and 
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mice versa. For the model developed Pee smetesi ss Only 


mbass one effector impairments are considered. 


mee CONCLUSION 

In this chapter the mathematical models for each compo- 
Meat in the block diagram of Fig. 2.3 were developed. fThe 
Me vidual models were then assembled to form the complete 
model of the F/A-18 dynamical system. Next the simulation 
program is developed to compute the model matricies in Eqs. 
3.54 and 3.55, and compute the response of the system to 


Peeeck and rudder inputs. 


> 


IV. PROGRAM DEVELOPMENT AND MODEL VARIA ea 


A. INTRODUCTION 

To validate the F/A-18 system model a computer program 
was written to compose the model matricies in Eqs. 3.54 and 
3.55, and compute the system response to Stick and rua. 
inputs. In addition the procran simulates an actuation loss 
of the right or left stabilator. The program was written in 
VS Fortran on the IBM 3033 computer at the Naval Postgraduate 
school. The program is organized to offer flexibility 


future development and modification. 


Ben PROGRAM, SiRUGiURe 
The program may be divided into four major operations: 

iy Dara aout 

2) thay idekearec ompiit aeaegis 

3) Composition of the system matricies 

4) System response computation 
The operations are performed by a series of subroutines as 
Shown in the flow diagram in Fig. 4.1. In Fig ee 
solid lines indicate control flow, and the dotted lines Sie 
cate data flow. The subroutines perform the steps outlined 
in Chapter II to compose the system matricies and compute 
the response. A £unctional description of each SuUbroimmmeme 
shown in Fig. 4.1 and associated variables is given in Appen- 
Gd tegen dS eee subroutines at NPS were used to perform the 
required matrix manipulations [Ref. 5]. These subroutines 


are also defined in Append: <a 
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mead in a/c stab. 
ema cont. 


Peabo DATA i 
ees. .f Light Alo 

eSnaitrans control and compose the total 
stab. & cont. failure system matricies 
Blerivatives CaaaNe Germs 


Taeeraccive 


LAWS | 
compose the 
Seneca law 
matrcies 


| 

| FLITE1 
| read in s.s. 
| 

| 


eOnaitions -& 
Sole eae ad Luce 
Perameters 


je ay 


compute the 
Variable gain 
vel ti 0 


Cem ce 


\ 


AIRDAT 


Pee ee 


compute discrete 
system equations 


sompute air datal | 
Een control law 
Pou 


ARCRFET 
compose airframe + 
actuator + sensor 
matricies 


| 
LONLAT | 


compose LONG & 
meri) mactricies | 


Compose unimpaired 
Ereerrame matrix 


MODEO PLANT 


compose modified ae oe 
eri ramne Matricies 


compose airframe 
a Gee Dla eene 
matricies 


| Sill col 


OU 


ae ne compose sensor 
sompose actuator matricies 


Mmacricies 


7 ees Daca Flow 


=" COnt@rol EF Low 


Figure 4.1 Simulation Program Flow Diagram 


Shey) 


Execution of the program is controlled by the” Pao ia 
program. The Exec jeenee the input and output data files as 
shown in Table 4.1. le weer options are provided to define 
data input/output files. The result of each subroutine com- 
putation is written to the FAILS RESULT file, sna eeaa 
viewed by the user after program execution. A copy of the 


Exec program is given in Appenci ee 


TAREE ee 


INPUT/OUTPUT FILE DEFINITIONS 


Record Device File Type Use 
Number Name 
Ol Disk FA18 DATA Contains the FA18 program 


input data 


O72 Disk FA1LS8 RESUESD Contains the rest 1, 


the operations performed 


by each subroutine 


03 Disk OP AMA DATA Contains FAI18 system 


Matricies (Eqs. 3.34.0 


3.55) formated forme 
control system design 
package at NPS. 


O4 Disk OPT ECR DATA Contains the system 
response data computed 
by the FA18 program 
formated for the inter- 


active plotting progam 
at: NPS 


1. Datei I rie 
Subroutines ‘FLITE1' and ‘FLITE2° performepiteme ee 


input operations. All data is read fromthe FAIS data sia 
except the control and failure parameters which are jimpan 


interactively. 


ie 


ST leemceddaonme the aircraft steady state flight 
Meomditions and control and failure parameters. As indicated 
[oer ie. 4,1 the flight conditions are read from the FA18 
MetA file. The control and failure parameters are read in 
[/imeeractively. This allows the user to conveniently run the 
Meogram for various control Pepuesceand comenrolesuriace fail-~- 
Mees. FLITE2' reads in the basic airframe stability and 
@entrol derivatives from the FA18 DATA file. The stability 
and control derivatives are arranged in the matrix format 
Meow in bas, 3.33 ~ 3.36, and in Appendix E. 

2. Air Data Computations 

Pun etomNGompesimo Lhe System matricies, air data 
computations are performed by subroutine '‘'AIRDAT'. The sub- 
memedne computes the air data inputs to the control law func- 
tions using the standard atmosphere equations [Ref. 6]. The 
program does not include the logic to compute atmospheric 
conditions above the gradient (troposphere) region. There- 
fore the computations are valid only up to 36,000 feet. 

53. Composing the System -Matricies 

As seen in Fig. 4.1 the system matricies are com- 
posed by ten subroutine operations. The operations and 
associated equations in Chapter [III are as follows: 

ie Compose the LONG and LATD matricies Eqs. 3.38 & 3.39 


Pee Compose the modified airframe matricies Eqs. 3.40 & 
ce | 


Pee Compose the actuator matricies Eqs. 3.42 & 3.43 


Sie) 


4) 


5) 


6 ) 


7) 
8 ) 
9) 


10) 


Compose the sensor matricies Eqs. 3.46 & 3.4/7 


Compose the airframe plus actuator Matricie oa. 
3.44 & 3.45 


Compose the airframe plus actuator sensor matricies 
Eqs. o948° 8° 3449 


Compute the discrete system matricies Eqs. 3.50 &@ 3.51 
Compose the GAIN matrix siqs reo 
Compose the control Law matricies Eqs. 3529 5Gaee 


Compose the total system maricies Eqs. 3.254 355 


The basic component matricies are composed by first generat- 


ing a null matrix of proper dimensions, and then assigning 


the coefficient values to the proper elements. These matri- 


eles 2nelude:, 


) 
2) 
3) 
4) 


») 


LONG and LATD 


Actuator 


Seneo rs 


Controle law 


GAIN 


The coefficient values to actuator and sensor model ame 


written into the program as constants; they are not contained 


on a separate data file. The coefficients in the LONG, LATD, 


control law, and GAIN matricies are first computed, epg 


assigned to the appropriate matrix element. 


tine 


'LONLAT’. 


The LONG and LATD matricies are composed by subrou- 


reflect control surface damage for one Of the four seomeem 


Surface impairment classes. In the present program only 
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As explained in Section I11.C these matiamegmee 


class one effector impairments are simulated. For EIC=1 the 
appropriate control surface coefficients in the 'LONG' and 
mpAtD”’ matricies are set to zero. 'LONLAT' contains the logic 
Beeimpair the right or left stabilator. ‘LONLAT' also con- 
poses the unimpaired airframe control matrix, GmO. This ma- 
Meas will be used to compute the impaired gain matrix in 
Smproutine 'VGAIN'. 
The control law matricies are composed by subroutine 

mews. The subroutine may be divided into four operations: 

ie Compute function gains 

e Compttte the filter coefficients 

eee Compute the matrix coefficients 


4) Assign coefficients to the control law matricies as 
imme as yikes DoZ. 


The gains are computed according to the function definitions 
Piven in Appendix A. The functions use the air data computed 
Mame AU RDAT', steady state AOA, and normal acceleration to 
compute the gains. All functions are programmed exactly as 
Shown in Appendix A except for the following non-linear 


mmcetkions: 


Piimeea on: 92.0 Pome tewdimale stick gradient 
Bine © em sl Lateral stick gradient 

is) phe 1c abteyigieeg Rudder pedal gradient 

Hine Gomes 2 RSRI non-linear gradient 


These functions were programmed neime Bien ene at mie mic tof a 
[eeror Series Expansion about the origin. Some of the func- 
tions computed in the program are not used in the control law 
model. They are provided for user information on system per- 
formance. These functions are: 


ull 


Funetion 24 Trailing edge flap schedule 
Fume tions Trailing edge flap schedule qc limit 


SS ONoMehe abe jel 2 7 Leading edge flap schedule 

Funetien 29 Leading edge flap schedule ri limit 
FF Ueno Mi eo Nz Limit on AOA feedback 

Is(iigkeieakeigt “Il Rolling surface limit schedule 


Function ll2 Lateral acceleration an 
Function 113 Lateral acceleration gain 


The filter coefficients were computed using the procedures 
outlined in Appendix B to transform the analos filterc uae 
digital filters. The control law coefficients were computed 
using the procedures Sete in Appendix C for transfoPmagme 
the control path transfer functions into state variable Goma 
The function and coefficient values are output to the Fae 
RESUEE pales 

Subroutine 'VGAIN' composes the GAIN matrix based on 
control effector impairment. For the unimpaired system the 
GAIN matrix appears as in Eq. 3.31. The subroutine ise 
signed to implement the reconfiguration algorithm described 
in Reference l. 

Prior to composing the total system matricres ane 
discrete form of the continuous system matricies must be 
computed. This is done by subroutine “EXPINIT” | Rete 


"EXPINT' computes the matrix exponential, 


Fos*ts 
a3 p 


and thegant ecr ak, 
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The discrete system equations arethen computed as: 


ee a ps*ts 
and 
ts % 
Bps =e PS "ds x Gps 


The Aps and Bps matricies are written to the OPTMATD DATA 
jmebe. The data file is formatted for the control system 


design package at NPS. 


dies Response Computations 


SV oremeresponse 1s Computed for the recursive equa- 
iO 1 


Bsr al Slap ee. Ss sk pepo ax L 
Gia tS A Es Melee bE Is Uk) 


[ime response is computed for 500 data points. The AF18 and 
BF18 matricies are defined in Eq. 3.54. The state and input 


vectors are defined as: 
X(h)" = [ Xm(k) | Xa(k) | Xs(k) | Xe(k) J] 
C 
U(k)° = [ px py pz ] 


Mespomsemdatantopeali 55 states is written into the 
PeIPLOT DATA file which is formatted for the interactive 
memeerone rowtines at NPS. The variables contained in the 


state vector which are relevant to the user for viewing are: 
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el = ef = x2] = edie 
<2 <6 tae pal XZ = Siew 
ee ania x oe slots x25 = dar 
x4 = thed pout ae P tebe x27 = cau 
OMe xl7 = dler x29 = Baie 
Oe = iP xl 9 =, died peo I ale IL 


The output equation (Eq. 3.55) was mot proceanmoar 


C,. PROGRAM TESTING AND MODEL VALID 
To test the simulation program, four Sets Of) runs 
made for the followine contre! snpmease 
1) Positive longitudinal stick, no failure 
2) Positive longitudinal stick, right stabilator {aan 
3) Positive lateral sticw. moe!) fad lures 
4) Positive rudder pedal, no failures 
All deflections were -!1 inch step inputs of 3 second duraueemee 
All runs were made at .6 mach/10000 feet. The sampling rate 
was set to 80 hz. (Additional runs were made Ae Ac hz sam- 
pling rates with no noticable difference in the output re=- 
sponse.) Aircraft steady state data, including the stability 
and control derivatives, are given in Appendix J. The iaame 
RESULTS file, including the function and control magn sia 
efficient values, is given in Appendix Kk. 
The model was verified for: 


1) Correct direction of motion of contro lest t ae] ere 
corresponding alrerart Mottram 


2) Expected aircraft response for a right stabilator 
failure. 


3) Proper augmented aircraft motion. 
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mererence / contains response plots for the aircraft for simi- 
Mmeyeetiaeht comditions. these plots were used to verify the 
model for the proper augmented motion. Some of the plots 
Were reproduced for the thesis and are given in Appendix L. 
The control inputs used in Reference 7/7 are slightly different 
than those used for thesis model. Some comparisons however 
ean still Bede All response plots for the thesis model 
aré given in Appendix M. 
Miemsiterart response EO a poSitive longitudinal input 

mea To failure conditian is shown in Figs. M-1 - M.3. 
ime motion of the right stabilator, and leading and trailing 
edge flaps compare favorably to the McDonnell model response 
meme in Fie. L.i. In the thesis model the right stabilator 
initially travels -2.5 degrees. After a .5 second transient 
period the stabilator continues to -3.6 degrees in 3 seconds. 
lime corresponding McDonnell response shows an initial travel 
of -2.2 degrees. After the .5 second transient period the 
Stabilator continues to -1.8 degrees in 3 seconds. The 
major difference between the two responses is the stabilator 
Beate , the thesis model being slightly greater than the 
Mewonnell model. (Possible explanations for the discrepancies 
in the model are given below.) 

The leading and trailing edge flaps, driven by angle of 
attacks feedback, deflect approximately 7.5 degrees in 3.4 
Seeenads as shown in Figs. M.J and M.2. The corresponding 


McDonnell response show the leading and trailing edge flaps 


oe 


to deflect approximately 4 degrees in 3.5 seconds. ~(Notemeae 
initial conditions of the flaps show Fig. Lol are nee 
as in Figs. M.1 and M.2.) Aircraft pitch rate andes 
angle response for the thesis model are shown in Figs. M.2 
and M.3. The thesis model achieved a maximum pitch rate of 
8 degrees/second compared to 4 degrees/second for the 
MODeeme model. Also the McDonnell model achieved a con- 
stant pitch rate in 2.5 seconds. The thesis model achieved 
a 2.5 degree pitch attitude in 5 seconds compared to 15 
degrees shown in.Fig. 1.1. The increased travel of the 
stabilator shown in Fig. M.1 (compare Fig. L.1) may expueue 
the discrepancies in the pitch and flap responses. 

Figure M.4 shows the response of the aircraft with the 
right stabilator failed. Note the decrease in piteh ame 
and pitch angle magnitudes. Figure M.4 shows a negative 
roll rate which is to be expected for a right) stapmileumem 
failure. Reference 7/7 did not give aircraft response for 
CONttol Sureace tanmiyres: 

Aircraft response to a positive lateral stick inpieeee 
Shown in Figs. M.5 - M.8. The motion Of tite stapaetapeee 
aileron, and trailing edge flap compare favorably to the 
McDonnell response shown in Fig. L.2. The aileron in the 
thesis model achieves a steady state deflection of 1.1 
degrees in .5 seconds. The McDonnell response shows a steady 
state deflection of 1.6 degrees in .5 seconds. The )t ieee 


model stabilator (Fig. M.5) deflects 0.45 deemeca ieee 
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seconds, the aire spe Teale McDonnell response shows a 0.6 
Meeree deflectiongin 0.5 seconds. The trailing edge flap 
in Fig. M.6 deflects 0.4 degrees in 0.3 seconds which com- 
pares favorably to the response shown in Fig. L.2. (Leading 
edge differential flaps are not used at the Mime ondivt 1 ons 
Peete hunction 93, which sets the gain on the leading edge 
map path, is computed as zero.) The rudder response Chae. 
M.6) shows an initial negative deflection which is expected 
for coordinating a right turn. The magnitude of the deflec- 
tions for the two models is approximately the same, however 
the thesis response is slightly more oscillatory than the 
McDonnell response. The aircraft roll rate, yaw rate, and 
fereeanigle response are shown in Figs. M./ and M.8. Again 
only slight discrepancies exist between the two models. 

ie aireraft response to a positive rudder pedal input 
Memsnownh in Figs. M.9 and M.10. The corresponding McDonnell 
response is shown in Fig. L.3. Comparison of the two models 
shows the shape of the responses to be approximately the 
same. However the magnitudes in the thesis model are very 
much less than the McDonnell model. For example the rudder 
achieves a maximum deflection of 0.07 degrees as shown in 
Meee M.9. Compare this to a 1.5 degree deflection for the 
meeonmell model shown in Fig. L.3. The response of the 
differential aileron and stabilator for the thesis model 
Memes. M.9 and M.10) is essentially zero. This is expected 


meee the rudder to rolling surface interconnect gain, 


oy 


function 39, is computed as zero for low angle =o) sapere 


(See Appendix A, function 39.) It is unknown why the aileroms 


and stabilators respond as shown in Fig. L.3. PossiDivywiea. 
response shown in Fig. L.3 is based on a different flight 
control program than that used in the thesis model. 
Possible explanations for the discrepancies (n@ gee 
sponse between the two aircraft models are: 
1) Improper derivation of the simplified contro. 
model based on the assumptions and procedures given 


Lin ap ber Sei a ee 


2) Programming discrepancies resulting in erroneous 
COnlputations: 


3) Errors in the aircraft linear small penprugeanwem 
model. 


4) Differences between the modeling techniques sa es 
in Reference 7 and in this thesis. (The MCGADK@iegem 
included all aerodynamic and control system non- 
linearities, as well as the effects due to digital 
time delay and quantization.) 

5) The McDonnell model response made available to the 
author may be based on a slightly different tlignhe 
program than used in the thesis model. 

The author investigated each of the items listed above. It 
is felt the assumptions used in developing the control law 
model are correct given the available information on the 
F/A-18 control system. The simplifying assumptions used vam 
developing the model may account for some of the differences. 
The program code was thoroughly reviewed and revealed no dis- 
crepancies. It is assumed however that this item remains a 


possible source of error. Two possible sources of error 


exist with respect to the aircraft péerturbdatvonmme a oe 
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fee ime coefficients in the NATC small perturbation 
model are computed with respect to the aircraft 
body axes (see Appendix E). This may effect the 
magnitude of the motion variables which are fed 
iaciememthneecontrol laws. 


2) The aircraft model developed in Reference 7 con- 
Siders the offset position of the accelerometers 
PaOlmmNeGmalrcratt center Of gravity. Lt has not 
been determined if this effect is considered in 
the NATC perturbation model. 


Further investigation into these Gees ba liters should be 


conducted. 
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IV. CONCLUSIONS AND RECOMMENDA 


A mathematical model was developed which simulates the 
dynamic response of the F/A-18 aircraft, Slite mogoueeee 
designed to implement the reciegie teams gain matrix fae 
the study of reconfigurable control systems. A program 
waS written which composes the discrete time state variable 
matrices of the alrerare. and computes one response to 
Strek and rudder inputs.- Die ppeenan also simulates the 
actuation loss of the right and left stabilators. Response 
plots of the thesis model were displayed and compared to 
the aircraft model developed in Reference 7. Possible 
sources of error were discussed, and it is recommended that 
further investigation into each of these areas be conduememm 

It is also recommended that the NATC program which com- 
putes the stability and control derivatives be acquaredwire 
the Aeronautical Engineering Department at the Naval Post- 
graduate School, to be made an integral part of this program. 


This would expand the ability of the program to simulate the 
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FUNCTION 24 


TRAILING EDGE FLAP SCHEDULE 
AUTO FLAP UP (AOA, Rj) 


(M s 0.61 Ry s 0.27 


20 
M2 st. iienti3.5. 17) 
IM = 0.9! Ry = 0.66 
SCHEDULED 
POSITION | | 
ae 16.4, 9) ia 3) 
(iM 2 1.05) Ry = 0.9! 
ine “See e. \125.6, 0} 
iG 10 20 30 
ALPHA - DEG 
a F24 
{(Oes} 


PEPriA 
(DEG) 


a ode 


CF} 


[428769 EAS 


Ce F24L1 


$2.76 12 AS 


63 


FUNC tf i UiN es 
TRAILING EDGE FLAP SCHEDULE 
AUTO FLAP UP (Qc LIMIT) 


ZO 


(6007. Oe ee 


ee ee re a 
O 512) 6 1000 L=U0 20Ge 


47, 63695145 


LIMi ft 
B25 


(eS) 500 


84 


ee Oi 2 / 


LEADING EDGE FLAP SCHEDULE 
AUTO FLAP UP (AOA, RT]! 


/ 4 
/ / 


Lil ensio) es 


40 FUNCTIONS 
28 AND 29U 


—- 4 7 
7 / 


(QaWS 2 == lee vommecrs 


G16; 


Sere OVRep 
pUSlPICN 20 
DEG 


0.0) 
ee one 10 20 30 40 
ee ALPHA - DEG 
LIMITED BY 


FUNCTION 2S 


ae ALB (DEG) 
Peese4) Bl AS 


85 


me Gr SIN 2s 


LEADING EDGE FRAP Sep eeem= 
AUTO FLAP US eee ae 


40 
[SO noo! (Zils Oger eo: 
30 
LIMIT 
DEG = 
1O 
(O50, een 
5 200 400 600 BOO |O00 


Qc - LB/FTS 


4422 Snes 


Cit 
S110) 
rAa\@ 


S| 
(BeGH 


86 


FUNG GIN Zs 


LEADING EDGE FLAP SCHEDULE 
AUTO FLAP UP (Ry LIMIT) 


40 


fo (UPPER LIMIT) 


CQIMAcH =sOnea1 10.7, 341 


DEG <0 
Lo 
Re 
ge) 
F2SL (LOWER LIMIT) -STORES OR STORES AND ais Se PGB RS carn ol 
/ [MACH eo bi.) 45, ©.0) 
r ; : F IN 2 Fd vo TD EAD ES ER A 
wu a ons eo | Were 2.0 
Ry UG=7Fe} 


eo ee OW ia Ml ty -Siores ANS a Zeal Zoosesh 


HP cae! = Wieser el fra 


F29L ISTORES AND Ps 2 1200 PSF) 


Ps 
(PSF ) 
F2eoL 
{DEG) 
Ry 
{Cc/Ps) 


O-OrSa AS 


F29L = O DEG FOR STORES OR STORES AND Poeun 1200, PSE 


87 


FUNC TO Nie 


DIFFERENTIAL TRAILING EDGE FLAP GAIN SCHEDULE 
AUTO FLAP UP {Ry}, Ps): 


O.2F 
(dl87 400-16) p= AOST eons 
(1.058, 0.18) 
10.518, 0.18) \ 
\ ee 
« e « ] 
1 ee (1.378, 0.12491 
GAIN a 

{VOLTS/VOLT) Po « S73 PSF 

10.717, 0.Ce64) (1.378, 0.0870) 
ale © 1485 po 

(1.379, 0.CS55) 
Pe > 2118 PS 

().37¢9. O2olcs) 
06 0.5 1.0 a 2.0 
Ry 
0.2247 BIAS 
ua Liwiy ; LIMIT CN 
(PSF) ee ee (VOLTS/VOLT | 
( ) 
0.737 BIAS 
es 0.112 S1AS as ae 
al 
Ry Ox) LIMIT 
a : 0.682 
Pe 


88 


FUNCTION 32A 


LONGITUDINAL FORWARD LOOP GAIN SCHEDULE 
AUTO FLAP UP (Qc) 


UZOGy 0.2) 
(241.7 KEAS} 


GAIN 


PASS 18) Ve sy 
SSS) | Gy 0S) | 


e S500 1000 =are) 2000 2500 
Oc - PSF 
LIMIT 
ee = 
ZO00 ~ ee F3CA 
pP 
(PSF ) 500 QKF 


89 


PUNC | Pines 


LONGITUDINAL FORWARG EQOr GAU Ne Sse. eaee—e 
AUT OF sae a OTHETAQEL se 


Pe Se soe 


10.9135, ©.59)! 


GAIN ooo 


(6e0. 6239710 Ps = !000 
(0.6043, 0.05) (0-790 au Oma 
(0.5398, 0.0515 
[6.02 40.2eom8 Ps 2 2000 


(O20, O48, 1590) CO. 7 2s 


(044729 ;20.05) (0. 7O2e0n Gr eo) 


70.023 Gia> 


NO, 


FF UINIE TON Ss: 
DIFFERENTIAL TRAILING EDGE FLAP 
GAIN SCHEDULE 
AUTO FLAP UP (AOA) 


O 
(D 


O 
O) 


O 
aN 


GAIN 
WO rs7 Ver i 


Of ———— |e TORS) Ceo. 2) 


Cree ie 
EVGETS7 VOUT | 


OL 


ge NCIS ONS ss 


LATERAL FORWARD LOOP Galil sep=2 sas 
AUTO FLAP TURP 3. e238 


GAIN 


VOLTS/VOLT 0-2 
7 ese) 
O15 0 (0 20 30 
Ney. = eae 
evsse HIAs 
Beeces A 
> 6 DEG 


AOA 0,5 BAS 


GAIN 
(DEG) rn ee | 2.0 fe {(VOLTS/Volem 
5 


oe 


Bene ON: 26 


AILERON GAIN SCHEDULE 
AUTO FLAP UP (Qc, Ps, Ry) 


Oe (REVERSAL POINT) SCHEDULE 
1t+.9, '!9e@l) 


20S0 
{2.S, !99@1) 


HINGE MOMENT GAIN SCHEDULE (Ry, Po! 


(0.42, (0.88, 
1.01 1.0} 


ae 400 PSF 


{0O.8, 0.682) aS y Ste 


\ 


EVOLTS/VOLT} 10.8, 0.45) 


Fie tizgits PSF 


REVERSAL GAIN SCHEDULE (Oc, Ca! 


O 0.5 1.0 
Ry 60 1Gg-200, 1.01 
pee GR 
page ivaTs~woarT) °:3 
Ry Limit 
: 0.8 


0.79950) A5 ——- 


-20C 
ines x 1c? oR an ae 
440 BIAS [.0 i 
ae iiag , PuPPER | ou | PuPPER | ou | 
= . i an Feats aie GAIN 
poe xX oceanic Soe oe oie Pep 
440 Vas |: bo | | sect | IVOLTS/VOLT) 


95S BIAS Ge IPSF I 


-260 [6.005 | L2:0Cs } 
Vfer7Bl B38/V6E.0-6 


S18 


BOINGION Ss / 


NZ LIMi qd GN AOA P22 hine. 
AUTO FEAR Us ita.) 


he 
EAS, (336) eeee 
GAIN 
(YCGRS7) Clef) 
Oi 


(SO eee. 


NZA 2G 


2.o.8 12S 


| 
SA s.0 Se 
(G'S) a 6 a (VOLTS/VOLT) 


94 


EUNGT TON 33 


RSRI GAIN SCHEDULE 
AUTO FLAP LP (AQAA, R]}, Ps! 


0.3 


(-$.0, 0.2071 


“tO Pe) 


“0.3 tto, -0.4¢5) x 1:8, -0.54! 
aay T aS ,$ 0.25 
118, -°0. 83) 
ae eee -0.991 


Ar 2 0.30 


e 


PUB, +t. ht 120.4, °5.2) 


er ean ee eeanennwnoeteatn2ee2e2e2222 


~- 138448 FSAI - .01774 


LPeER 
Limi ft 
120.4 
Limit 
ADA BAIN 
10€0) ay —O— WO TSvoLT) 
0.048 BIAS ADA ¢ 10 DEO 
0.0498 7 
Limit Limit 
ra ° ‘ 
+300 be} 0.01 beY>® 2.5 . me 
lat £240 0.0 y - 
0.72 BIAS 3.9 BIAS 2.3 BIAS 
Liwit 
mn 62) “210s --[orores 
0.0 


Js 


FUNCTION 39 


RUDDER PEDAL TO ROLL CAS INTERCONNECT 
GAIN SCHEDULE 
AUTO FLAP UP [AOA) 


ce) 


GAIN 
VOETSZ 7 eli: 


IMIT 
AOA hs ease eS GAIN 
(BEC i Peete i iVOLTS/7 Velen 


ee 


96 


FUNCTION 40 


PlienerhAve f SECO AEX GAIN SCHEGULE 


AUTO FLAP LP (Ps, Oc, RR]! 
0.8 


10.0, 0.&38) 
0.@ 


10.0, 0.475) 


| ae 


Bee 9721 een, 6.21) 
: 6.2 7 i 
: 11900, 0.1! 
11900, 0.07) 
0 2500 
Po > Por 
3.235 BIAS 
A, < 0.75 
ny Ay 2 0-75 
0. @58625 BIAS 
Limit 1 wo R 
ae °0.20177 BIAS o 
$00 
9.635 x 10°* et 
Py 
(rer) 
0.473 BIAS 
LimtT iprg on 
1 BOO -@.5 x 10°% x 
0.0 — 
1.5238 x 10°2 Bias ° PIQ ¢ S00 
“ Axe 
235 BIAS ar os 
Limit |~ 
es 108 
ipo) 
0.0 S$.0746 x 10°* BIAS 
1.9m) x 10°7 5 ives 
ae Q- < 440 


ome 


OX 


1.87247? x 10°? Bas 


: > 


oy 


AOA 
{DEG) 


Ry 


MENG EIN G | 


ROLLING SURFACE LIMIT SCHEDULE 
AUTO FLAP UP 1AOA,R]) 


LIMIT 
(VOLTS) 
O. 
@.0 BIAS 1.0 BIAS 
LIMIT 
+22 0.052571 eS 
2 wae | SS VOUT! 
1 
LIMIT 
+0. 48 0.23533 CX) bes 
tO oa? ASA g 8.0 
- 009387 
0.5875 
UPPER 
rt -—[s}- 
0.39 as 


he 


me NC) ONT 42 
RSRI NONLINEAR GRADIENT 


" 
se 


Pease 5U)) 


Ze 


SUTeUT 
VGLiS 


X10.0479/x/+0.7611 


, x 
- 30 10 20 
INPUT - VOLTS 
-10 
-20 
(-18.33, -30)} 
-30 
O.76! BIAS 
INPUT eae ABSOLUTE Palen | Se OUTPUT 
(VOLTS) $18.33 VALUE ae \ Poles) 


oo 


PONS] | ON 45 


DIRECTIONAL FORWARD LOOP GAIN SCHEDULE 
AUTO FLAP UP (Qc, Ps! 


ES § 625 PSe 


GAIN ee 2 1450 PSF 
(YOLTS7VOET} ; (1000, 0.5) 


0 2a0) 5 


750 [O00 i258 


| -O Vets 


7 oy GAIN 
/ N/" (VOLTS/VOLTI 


100 


FUNCTION 6&8 


PITCH RATE FEEDBACK GAIN SCHEDULE 
AUTO FLAP UP (Qc) 


E45 On ©sosa) 
Ca eas? 


mieis| ORG, 


o2StOl 5 eR, 
[S60 Rn GrS' 


"9 100 200 300 400 $00 
Qc - PSF 
-480 BIAS 


F68 


-2.9773 xX 1073 


JO) 


FUNC TONS Se 


LATERAL ACCELERATION FEEDBACK 
GAIN SCHEDULE (Ry, Ps) 


40 
l(2.8,> 35-0) 
30 
pe 
aia (2.8, 24.75) 
PS 
_ 435 Eel 
GAIN 20 ine 
VOLTS/VOLT Pe > 628 PSF 
(0,725 1S. 5 
lO 
O 
O oO. 5 ©) bes Fg © a5 


GAIN 
VOL IS7Vear 


O27 2 Bins 


Oe: 


FUNCTION S3 
OIFFERENTIAL LEADING EDGE FLAP GAIN SCHEDULE 
ACT CCmeArP Ur RT, PS]» SUL A= | 
(0.58, 0.06) (0.85, 0.06) 
Slots 10.434, 0.08} 10.74, 0.06} 
0.04 ee 
S| 8] 8} 8 
FORA ny a 
0.02 a Ae ae Ae 
10.42, O} (O-51, 0) 
[0.294, O) a 10.5 0} 
0 Gee 0.4 0.6 0.8 Rs 2 
= R; 
5 I7@8S, {1.01 {1000, 1.0) 1 gelO-Oe 1-0) 
F938 0.5 Fe3c 0.5 ines hore ee 
0 
0 8.0 
Ry 
Ps] 
(PSF } 
628 BIAS 
3 GAIN 
SQ is ete 
We Suh F93C 
(FROM FIGURE 11 es cae ~s 
1.0 BIAS 


ee 


BONE | LON - POC 


SPEEDBRAKE COMPENSATION INCREMENT 
AUTO FLAP UP (R}7) 


HO Sy Ziel, (0-85, 2.08 


BELA 
(eat) 


O @s2 OR, 0.56 O.8 ‘wae le 


Ry $ 0.7733 


19.2 BIAS Fr toe 


Ry > O, 72783 
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PIGING 1 CIN] Sey 


Ciena Nee Sen ene nen LEYIAT!ON SCHEDULE 
AUTO FLAP UP IR], Ps, NZ,e) 


Nz = 106 10.91, 0.015) 


GAIN pee 
tVOLTS/VOLT! ~° 10.78, 0.0091 
Ps 2 1339 PSF By ne 
Fo. = (OGG sr 
10.58, 0.01 
0 
: - uo 
Ry 
0.58 BIAS 
LOWER 
Ry Ox LIMIT 
a2 0.0 
Limit UPPER 
2 oe was 
MIT : 
Ps (ox | +835 LIM] 5.4446 xX 10 EX Pe 


= \o! 


498 BIAS Lio! 


798.54 5.A5 


Ws 


FUNGTIeNTSs 


YAW RATE GAIN SCHEDULE 
AUTO FLAP UP (Qc, Ps} 


(y34.1.1) Pe S “42 1670027) 


°o 200 400 600 BOO Kovors) 1200 
Qe - PSF 
1.54 BIAS 
UPPER 
LIMIT 0.0044 CX) 
: 168 Qe <¢ 500 PSF 
LPPER 
ae 0.38 BIAS S655) Hii GAIN 
(PS ; é ae (VOLTS/VOLT I 
Qe > SOO POF 
0.00089 CX) 
628 BIAS 
Pe LIMIT I 
+828 - 0 .COO77 72 
(PSF) ese as 
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eG ONP | O 7/ 


LONGITUDINAL INERTIAL GAIN SCHEDULE 
AUTO FLAP UP (Qc) 


GAIN 
(VOLTS/VOLT) 
ee] 
Og WOO B00 1200 1600 2000 2400 
Qe - PSF 
750 BIAS 8.4 x 10°4* BIAS 
GAIN 
(VOLTS/VOLT) 


WG) a 


FUNC?) TONS es 


DIRECTIONAL INERTIAL GAIN SCHEDULE 
AUTO FLAP UP (Qc) 


0.016 


(50770. oles) 
0.0le2 


GAIN 


(VOLTS/VOLT } ee 


0.004 


[39S0, O60 )e=) 
E200, 0. Ceres) 


12Q0 


GAIN 
(VOLTS/VGE 
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PUN | Beni 


RUDDER PEDAL COMMAND GAIN INCREMENT 
AUTO FLAP UP [R]) 


GAIN 
(VOLTS/VOLT} 


[e5o, 0.0) 


GAIN 


R 
(VOLTS/VOLT} 


0.14286 


SORE: 


FUNGIgIONe ae 


LATERAL ACCELERAT TEN eAty 
AUTO FLAP UP Ags 


1.00 (18.00,1.00! 


GAIN : 
(VOLTS/VGLT! 


(12.00,0.00) 


-00 4.00 8.00 t2.00 18.00 20.00 24.00 


AGA HEC 


AQA 
(DEG) 


©. e657 


GAIN 
(VOLTS/VGaw 


Mee ENT x Bb 


PRGITALD FILTER MODELS 


fee PITCH RATE LEAD-LAG FILTER P2 


.OI5(1+F22)S + 1 
P2(S) = ------------------ (B.]) 
SOUS S: 4 1 


Dsimg the tustin transform let 


2 (7-1) 
a 
ts(Z+1) 
which gives 
(P2N1 ) (P2N2) 
2{,. O15) (14+r27) + ts Somme c te Ola) (1tr22) 
een a 7G ee HK -————— | 
2(.015) + ts Lies) oes 


2{.015 )-+-ts 
(B.2) 
(P2D) 


RB. NORMAL ACCELERATION LAG FILTER P5 


P5($) = -----~---- (Bhs 3h 


Using the Tustin trans tome 


(PSN1) (P5Ne2> 
ts 
Se eee i 
O8 + ts O8 + ts 
P54 Zi) = «62256222 
O08 - 
oe || SSS aoe 
O8 + ts 
(PSD) 
Ce INTEGRATOR P92 
i 
P9(S) = --- 
iS 


(B.4) 


(PoSF 


Reference 3 approximated this filter using the backward 


difference method. Let 


Z-1 
S = ----- 
tsZ 
which gives 
(PON1) (PON2) 
ts7 + 0.0 
P9(7,) = ------3-----C-O 
Z - ) 

(pon) 


le 


(R.A) 


fee LEADING EDGE AOA LAG FILTER Pll 


] 
Pll1(S) = ---------- 
Oes0S + 1 
Meeng the Tustin transform: 
(P1l1IN1 ) (PLIN2) 
Es as 
| SO (4 eee 
mote) + ts .39(2) + ts 


.39(2) + ts 


(PLID) 


meet rALOING EDGE FLAP AOA LAG FILTFR P12 


Using the tustin transform: 


(P12N1 ) (P12N2) 
ts ts 
(ee 5- > se Yash, =o la ec 
~79(2) +° ts P/O ee ts 
P12(Z) = ---------------------------- 
70(2) - ts 
To ae cs a 


(P12N3) 


ies 


(Bee?) 


(B.9) 


et ©) 


(B.11) 


FF. YAW RATE CANCERLERV Ee 


(Y3N]1 ) (Y3N2) 
2 -2 
[------ JZ + w------ 
2 + “ts 2+ ts 
NAO) ee SS ae Se Se 
2-its 
Tu = Saaz. 
2+ ts 
(Y3D) 
G. RSRi LEAD LAG FiLtER Yo 
75S + ] 
¥S(S) 2] ——————— —— — 
5S + 1 


Using the Tustin transform: 


(YSN1] ) (YSN2) 

2( 75) 68s ts = Zeas. 

= = ee JZ + s-2-e------ 

2 25)) ees 2(.5)) +s 

YS(Z) = =-==2252 225 oe 
2\ 5 ees 
y eS, Me 
2(. 5 )eeeets 

(YSN3) 


114 


(B. 128 


(B. le 


(B.14) 


(B. ae 


APPENDIX C 


SIGNAL PATH TRANSFER FUNCTIONS AND STATE SPACE MODELS 


feePiTCH RATE TO COLLECTIVE STABILATOR PATH H1(Z) 


Peom Pra, 3.1 the forward path transfer function is 


H1(Z) = {P9*P2*F68*F32a*F12 + P2*((F68*F32a)+F40)1 (C.1) 


Let 
Re= sos fe 372n hi Z 
B = (F68*F32a)+F40 


Rewriting Eq. C.l using the filter Z-transform expressions 


for P9 and P2 gives 


A*(PON1*Z + PON2)*(P2N1*Z + P2N2) 


(Z-P9D)*(Z-P2DP) 


B*(P2N1*Z + P2N2) 
i eS 2 or ee ee iG) 


Bepanding Eq. C.2 


ed Z ) 


{fA*PON1*P2N1 + B*P2N1] Z 
+ [A*(PON1*P2N2+P9N2*P2N1]) + B*(P2N2-P2N1*P9OD)] Z 
+ f[A*PON2*P2N2 -— B*P2N2*PODI]} / (Z-P9D)*(Z-P2D) 


(Gas) 


eS 


Equation C.3 can be expressed in terms Of the fee. 2a 


partial fraction expansion: 


2 
bO*Z + bL*Z + b2 opsiell austen 2 
H1(Z) = ---------------- = qst3 + ------- + ------ (Caae 
(Z-P9OD)*(Z-P2D) (Z-P9D) (Z-P2D) 
Where 
ast3 = bO 
y, 
astl = (bO*P9OD + b1l*(P9D) + b2) / (P9N-P2D) 
2 
qst2 = (bO*P2D + bhl*(P2N) + 62) / (P2D-P9N) 


Equation C.4 can now he put in the following state variable 


forms: 
x1(k+1) POD O } xl(k) ] 
= A ne) (Cas 
x2(k+1) O P2DIix2(k) if 
rag le 
estxl(k) = [ast] ast 2] + ast3 atk (CoGe 
Rois) 


The sytem of nomenclature used for the coefficients in the 
output equation (Fq. C.6) 1s described imececemieee. 


In the remainder of the appendix only the first and last 
steps of the state variable development will be given for 
each transfer function. The equations for each of the 
coefficients are given in the simulation proaqram listina. 


ee) 


Nea AlwACeCELERATION TO COLLECTIVE STABILATOR PATH H2(Z) 


H2(Z) = ([F12*F32a*3.5*P9O*P2 + 3.5*F32a*P5 | (eae) 
x3(k+1) POD OO x3(k) 1 
= + nz(k) (c.8) 
x4(k+1) O P5D x4(k) ul 
a LK) 
Setx2(k) = nzstl nzst2 Homes ts oc i ) Ceres) 
x4(k) 


PITCH STICK TO COLLECTIVE STABILATOR PATH H3(Z) 
H3(Z) = [7.0*F32a*F12*P9 + F20*F32al Genie. 


Note: F20 defines the longitudinal stick gradient. The 
value 7.0 is a taylor series approximation 
about the origin of F20. (See F20 in App. A) 


x5(k+1) = POD x5(k) + 1 px(k) Cea) 


Pee) — OXStl ~5(k)) + pxst2 px{k) CCEsi 24 
ANGLE OF ATTACK TO COLLECTIVE LEADING EDGE FLAP H4(Z) 
Beg — 1.328*Pil reales.) 


Note: The leading edge flap schedule is defined by F27. 
The value used in Eq. C.13 is the l.e. flap 
gradient. (See F27 in App. A) 


x6 (k+1 ) 


Bie Dexouk) + 2 aalk) (ela) 
eikexi(k) = aalel x6(k) + aale2 aa(k) CGS) 


ele? 


ANGLE OF ATTACK TO COLLECTIVE TRAILING EDGE FRAP Hoe 


H5(Z) = 124054212 (Cc. lay 


Note: The trailing edge flap schedule is defined hy 
F24. The value used in Eq. C.16 is the t.e. flap 
gradient. (See F24 in App. A) 


x7(k+1) 


P12D x7(k) + iPaq (C. 13a 


eatex(k) = aatel x7(k) “=Seaate2eeae (Cc. ey 


ROLL RATE TO DIFFERENTIAL STABILATOR PATH H6(Z) 


H6(Z) = rv7*(F44+rk6t) (ci 
estyl(kj>= v7 (R42 kot) eerie (Cl2Gh 
where: 


rv7 = MIN(F6*F35,F6-F1]01) 
rk6t = LIMIT(O,.12-F4,.12-F4) 
LL UL Xx 
LATERAL STICK TO DIFFERENTIAL STABILATOR PATH A7tz 


H7(Z) =. 32224 rv7* 8/7 oe (Cc.2a8 


Note: The lateral stick aradient is defined hy Fl. 
The value 3.22 is a taylor series approximation 
ahout the origin of Fl. (See Fl in App. A) 


esty2(k) = 3.22*rv/*F7*{ P23 eae (Co22e) 
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RUDDER PEDAL TO DIFFERENTIAL STABILATOR PATH H8(Z) 


Bets )e= FEA“). 33*rv7*F39 Ger 3)) 


Soro Gees eet ose rv 7 FSS pak) fe. 24 ) 


ROL RATE TO DIEFEERENTIAL LEAPING EDGE FLAP PATH H9(Z) 


H9(Z) = F93*(F44+rk6T) er-25:) 


eleyl(k) = F9O3*(F4+rk6t) rr(k) Ve. 2G) 


PewERAL STICK TOeDIFEFERENTIAL LEADING EDGE FLAP H10(Z)} 


H10(Z) = 3.22*F93*F7*(F13+F4) (Cx 7) 


eley2(k) = 3.22*F93*F7*(F31+F4) py(k) Gezo8) 


ROLL RATE TO DIFFERENTIAL TRAILING EDGE FLAP PATH H11(Z) 


H11(Z) = F31*F34*(F4+rk6t) ceo) 


eteyl(k) = F31*F34*(F4+rk6t) rr(k) OUNEO)) 


LATERAL STICK TO DIFFERENTIAL LEADING F.DGE FLAP H12(2Z) 


Be wee ys — 93.2277 *P31*r34*( Fi 3+F4 ) (ery) 


etey2(k) = 3.22*F7*F31*F34*(F13+F4) py(k) (GE32) 


ik? 


ROLL RATE TO AILERON PATH H13(Z) 


rele eZ) P3957 F36 So Seo (C. 32m 


eal (k) F35*F36*.5*(F44+rk6t) rr(k) (C, 34a 


LATERAL STICK TO AILERON PATH H14(Z) 


H14(Z) 3.22*F7* F35*F36*.5*(FL3+F4) (C. 3a 


ea2(k) 3222727 Fos ow (Cl 3ee 
RUDDER PEDAL TO AILERON PATH H15(Z) 

HIS(Z ) = f14*Fs0* 1 .33"rs5 scG (C. 398 
ea3(k) = Fl4*F39*1.33*F35*F36 pz(k) (Cc. 3eH 

YAW RATE TO RUDDER PATH H16(Z) 

H16(Z) = F45*F96*COS(alpha) *Y3 (Cc. 3a8 


Note: In the simulation program alpha is set to the 
steady state angle of attack in deqrees. 


x8(k+]1) = Y3D xP(k) + 1 yr(k) 


erl(k) = yrril x8(kj)0+} yer2 ew 


L2e 


OQ. 


i. 


RPOUL RATE TO RUPERR FATA H17(Z) 


H17(Z) = F45*F96*Y3*SIN( alpha) 


+ F45*P38*Y5*F30*0.76¥*((2*rra)+(2*rrst)) (C.40) 


Note: The second term on the r.h.s. of the equation 
results from the rollino surface to rudder 
interconnect path. The value, 0.76, represents 
F42, the RSRI nonlinear gradient. It is a 
taylor series approximation about the origin 
of F42. 

x2(kK+1) 5 66))) e) x2(k) iL 
= + tri) (e741) 
x3(k+1) 0) Y5D x3(k) iL 
x2(k) 
See rrr =rrr2 erro rr ik) Cee4 2) 
oui) 


He (Z) 


LATERAL ACCELERATION TO RUDDER PATH H18(Z) 


BAST ho KG.430) 


er3(k) = nyr ny(k) (ome: 


S. LATERAL STICK TO RUDDER PATH = aa 


H19(Z) = F45*Y5*F38*F30*0.76*((2*pya (2 oy oe re 


x4(k+1) 


er4(k) = pyrl x4(k) + pyr2 py(k) (ee 


T. RUDDER PEDAL TO RUDDER PATH H20(7) 


H20(Z) = P45*Y5*F38*0.76*F30*((2*pza 2 ao ae 


- F45*F14*( .5-(F17*F114) ) (ea 


Note: The first term on the r.h.s. of the Gquatver 
results from the rudder to rolling surface 
interconnect path. 


x5(k+1) YSD x5(k)o + Sipteziia 


er5(k) = pzrl x5 U6) ets o2r2 area 


YSD x45) ee) Ge. 


.45) 


46) 


47) 


48) . 


x9(k+4+1) 
x LO(kK+1) 
xlL1(k+1) 


12(k+1)} 


[ x 


APPENDIX D 


CONTROL LAW MATRICIES 
[Ac] 
PON 
P2N 
POD 
PSN 
PID 
= Prep 
PlZzD 
we 
[Beet 
So ae: 
1 0 o! 
QO 1 Of 
Om 1. oO 0 Gah) 
e160 Yo! ‘ nz(k) 
OF OP aal(K) + 
Ce ee | 
TY “O07 OF lrr(k) 
Ch OL gies) 
0 | 2 1 O = 
e OO nO 
IV jee 1) 


Las. 


ron 


eS! 


[99] 


(%) Au eS 0 
9 ee (xj4al |o eaas- go | 2 
(>)2d O  {Ad| 0 y (HM) aA See eee () 0 i 
(CoO) WCayAd} faszd _asAd (4)eR ‘Bleek = 0) | 
as? = 0 (4) zu FG BZUL OS cun 
oO 7 0 4) b 0 CORR 0 
~ ysxd- | 0 ejSZU-¢4 sb 
[od] [95d] 

(4) 21x 

(4) TUX 7 
Seer iezdsyaAd gai4 4422 1124- (%)4a 
(4) 6x i) it a lk (ive 
(4) 8x | : (y)Aaqa 
(4) 2x x O (4 WvAGoe 
(4) 9x 6 a es es En ee 2 ES ee 
(Ug (Uae 6} 0 0 0 0 0 (4) xoqo 
(4) 7x O  1eFLLe O 0 0 0 O (4) xa,a 
(3) EX | O O LjASNG= %2 FS4uU 7145S 2U-7 955") Sb (%)xqSa 
(4) 7x 

( 
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APPENDIX E 


Seixas loiItTyY AND CONTROL DERIVATIVE DEFINITIONS AND UNITS 


The matricies in the small perturbation model (Eqs. 3.32 
through 3.36) are defined in this appendix along with the 
dimensions of the stability and control derivatives. This 
information was obtained from the flight systems hranch at 


the Naval Air Test Center. 


FX MATRIX 
l 2 3 
Ht <0) XW XO-WR 
74,0) ZW (.ZO+UB) 
ea 
1-ZWD vane 1-ZWD 
MWD* ZU MWD*ZW MWD* (ZO+UB) 
3 MU+------ MW+------ MQ+----------- 
1-ZWD 1-ZWD SA ago) 
4 O O iL 


ez 


A 


-G*CTHT 


1-ZWD 


MWD*ZD 


S) 


GX MATRIX 


HX MATRIX 


DX MATRIX 


1, 2 B 

O 0 O 
ZDS ZOE ZDTE 
1-2WD 1-ZWD 1-ZWD 

O O O 


JSea AS Nap eNah sc 10).4 


1 Zz 5 - 
YV YR-UB YP+WB Grea 
1-YVD 1-YVD 1-YVD }-YVD 
Ny Vay NVD* (YR-UB) NVD* (YP+WP ) WSC. enn T 
NV+------ NR+----------- NP+-----------  --------- 
1-YVD 1-YVD 1-YVD 1-YVD 
By D*YV LVD* (YR-UB) LVD* (YP+WB ) LVb*G*CnuT 
LV+------ LR+----------- LP+----------- 9 ---------- 
1-YVD 1-YVD 1-YVD 1-yYvP 
Sine 
0 ---- 1 0 
Cle sbk 


ea 


GYZ MATRIX 


1 2 3 4 5 
YDHT YDLP YDTF YPA YDR 
1-YVD 1-YVD 1-YVP ]=-YVD 1-YVP 
NVP*YDHT NVD*YDLF NVD*YDTF NVD*YDA NVD*YDR 
Vo NOLE.====—— Not? = 1 NDR+-=——— 
1-YVP 1-YVD 1-YVD 1-YVD 1-Vae 
LVD*YDHT LVD*YDLF LVD*YDTF LVD*YDA LVD*YDR 
LDHT+------- La pouor Se = — LPA+------ LD2+-— = 
1-YVP 1-<Yr 1-yvp 1L-YVP 1-YVP 
O 0 0 O 0 


HYZ MATRIX 


iE 2 5: 4a 
] O l O O 
2 O O {) O 

yo YR YP O 
; ies wo ae 


EAS 


DYZ MATRIX 
1 Zz 3 
O O 6) 
0 O O 
Yous YDIGE YDIF 
cr J—-YVD 1-YVD 


i =) 

O 0 

0 O 
YDA SDR 
1-YVD-1-YvD 


The following abreviations were used in the ahove 


matricies: 


BODY AXIS LONGTUDINAL WIND 


BODY AXIS VERTICAL WIND 


COS(THETA) 


SIN( THETA) 


DEFINITION OF D S ERIVATIVES 


R 


RIYATI 


LONGITUDINAL FORCE 
LONGITUDINAL FORCE 
LONGITUDINAL FORCE 


OF; 


LONGITUDINAL FORCE 


LONGITUDINAL FORCE 
LONGITUDINAL FORCE 


LONGITUDINAL FORCE 
LONGITUDINAL FORCE 
LONGITUDINAL FORCE 


VERTICAL FORCE 
VERTICAL FORCE 
VERTICAL FORCE 
VERTICAL FORCE 
VERTICAL FORCE 
VERTICAL FORCE 
VERTICAL FORCE 
VERTICAL FORCE 
VERTICAL FORCE 


PITCHING MOMENT 
PITCHING MOMENT 
PITCHING MOMENT 
PITCHING MOMENT 
PITCHING MOMENT 
PITCHING MOMENT 
PITCHING MOMENT 


PITCHING MOMENT 
PITCHING MOMENT 


LATERAL FORCE 
LATERAL FORCE 
LATERAL FORCE 
LATERAL FCRCE 
LATERAL FORCE 
LATERAL FORCE 
LATERAL FORCE 
LATERAL FORCE 
LATERAL FORCE 


ROLLING MOMENT 
ROLLING MOMENT 
ROLLING MOMENT 
ROLLING MOMENT 
ROLLING MOMENT 
ROLLING MOMENT 
ROLLING MOMENT 
ROLLING MOMENT 
ROLLING MOMENT 


es 0 


WITH RESPECT TO: 


FORWARD VELOCITY 
VERTICAL VELOCITY 
PITCH RATE 
VERTICAL ACCELERATION 
SPEED BRAKE 


THROTTLE 


HORIZONTAL STABILATOR 
LEADING EDGE FLAPS 
TRAILING EDG: FLAPS 


FORWARD VELOCITY 
VERTICAL VELOCITY 
PITCH RATE 

VERTICAL ACCELERATION 
SPEED BRAKE 


THROTTLE 


HORIZONTAL STABILIZER 
LEADING EDGE FLAPS 
TRAILING EDGE FLAPS 


FORWARD VELOCITY 
VERTICAL VELOCITY 


PITCH RATE 


VERTICAL ACCELERATION 
SPEED BRAKE 


THROTTLE 


HORIZONTAL STABILIZER 
LEADING EDGE FLAPS 
TRAILING EDGE FLAPS 


LATERAL VELOCITY 
LATERAL ACCELERATION 


YAW RATE 
ROLL RATE 
AILERON 
RUDDER 


DIFFERENTIAL LE FLAPS 
DIFFERENTIAL HORIZ STABS 
DIFFERENTIAL TE FLAPS 


LATERAL VELOCITY 
LATERAL ACCELERATION 


YAW RATE 
ROLL RATE 
AILERON 
RUDDER 


DIFFERENTIAL LE FLAPS 
DIFFERENTIAL HORIZ STABS ~ 
DIFFERENTIAL TE FLAPS 


UNITS 


yl FT/SEC2 


FT /SeG@e 

RAD/SEC2 
FT /Se@g 

RAD/SEC2 
PCT /SEG2 
RAD/SEC2 
RAD/SEC2 
RAD/StC2 


FT /SEGz 
FT /SEG2 
RAD/SEC2 
FT /SEGZ 
RAD/SEC2 
PCT/SEC2 
RAD/SEC2 
RAD/SEC2 
RAD/SEC2 


FT/SEC2 
FT /SE@a 
RAD/SEC2 
FT /SEGa 
RAD/SEC2 
PCT/SEez 
RAD/SEC2 
RAD/SEC2 
RAD/SEC2 


FT /SEGzZ 

F 1 /SEGz 

RAD/SEC2 
RAD/SEC2 
RAD/SEC2 
RAD/SEC2 
RAD/SEC2 
RAD/SEC2 
RAD/SEC2 


FT /SECZ 

FT/SEGZ 

RAD/SEC2 
RAD/SEC2 
RAD/SEC2 
RAD/SEC2 
RAD/SEC2 
RAD/SEC2 
RAD/SEC2 


NAME 


| OF DIMENSIONAL STABIL RIV 
(Continued) 


PARTIAL DERIVATIVE OF: 


YAWING MOMENT 
YAWING MOMENT 
YAWING MOMENT 
YAWING MOMENT 
YAWING MOMENT 
YAWING MOMENT 
YAWING MOMENT 
YAWING MOMENT 
YAWING MOMENT 


Iesyil 


WITH RESPECT TO: 


LATERAL VELOCITY 
LATERAL ACCELERATION 
YAW RATE 

ROLL RATE 

AILERON 

RUDDER 

DIFFERENTIAL LE FLAPS 
DIFFERENTIAL HORIZ STABS 
DIFFERENTIAL TE FLAPS 


ES 


UNITS 


micouea 

Ri obee 

RAD/SEC2 
RAD/SEC2 
RAD/SEC 2 
RAD/SEC2 
RAD/SEC2 
RAD/SEC2 
RAD/SEC2 


APPHENE IX 


ACTUATOR TRANSFER FUNCTIONS 


Al SSTABILATOR 


S 2 2(0.0629 


(----- ) + -------- See 
Dst(S) 82.9 82.9 
Est(S) S 2a Oa) SZ 2( 0; Sicm 
[ (----- ) + ------- S + Glee ) + ------- Sie 
36.4 36.4 LOS 1052.6 


The procedures outlined in Ref. 8 were used to put the 
transfer functions Im State varliableoweennle 


1) Arrange the transfer function as 


2 
b2 S + b3 S + ba 
A 3 2 
S +alS+a2S +a3S + a4 


D2 =“ 7s 7 E03 a2 = ]1.6122E+04 
b3 = 2.4101E+04 — a3 = 4.9559F+05 
b4 = 1.4691 E+0/7 a4 =81 746° OF 
al = 1.5410E+02 
2) Express the transfer function as a differential 


equation 
Dst(t) + al Dst(t) + a2 Dst(t) + a3 Dst(t) + a4 pea 


= b2 Est(t) + h3 Bsttt) => bau heeaee 
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3) Compute the state eauation coefficients 


UO = dO = 0 

tie —s ol —  al“~U0O =O 
Uze— 62 — al*Ul —"a2* 
We = 530-7 a!l*7U2 = a2x~ 
U4 = b4 - al*U3 - a2* 


UO 
OLE 
B2 


2.37 7E+O3 


a3 7~UG. = 


se ODS Z Ero 


a3*Bl - a4*BO = 2.7277E+07 


4) Arrange in the followina state variable form 


Fst 
peal ) i 0 
x2 = 9) 0 i 
ma @) 9) 0 
x4 -a4 -a3 -a2 -a 
Hst 
xl 
Deter 1. © 16 - «2 
x3 
x4 


xi 
xe 
x3 
x4 


Gst 


Ul 
+ U2 
Us 
U4 


BS t)) 


Note input and output for all actuator models is in 


degrees. The remaining state variable models are developed 


as above. 
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B. LBADING ELGE- FEA 


Dle(S) pe ©) 
Ele(S) S S 
(~--- + 1.0) ( ---- + 1.0) 
26.9 82.9 
a relia 
= Fle + Gle Ele(t) 
ae ic x2 Ce 
x Miey) 
Dieter) 2 erie 
Seo (ee.) 


GG.  ORATLEING EDGE EL 


Dte(S) noe | 
Ete(S) S 2 21 Oa) a 
(----) + ------- S + 1.0 | 
S50 35-20 | 
| 
ne) itt) 
= Fte + Gte Ete(t) 
eo Ce) x2UL) 
Sete) 
Dele) ee ee 
So 2 Gs) 
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Dee AILERON 


Da(S) liens @. 
Fa(S) = S 2 ~~ 2(0.59) 
(Ge ee eS S + 1.0 
Few 8, 75.20 
cole t ) Sole te) 
= Fa + Ga Ea(t) 
me ( €} x2(t) 
xt) 
Wart) = Ha 
ef es) 


Eee RUDDER 


Dr(S) LAG 
Er(S) S 2 2(0.69) 
(552) 95 = eeeee S +1.0 
Ook ee a 
late ad 6 } 
=| tte + Gr Er(t) 
x2(t) ote) 
cle te) 
Dig (te Sete 
Met) 


IL3 3 


The actuator state variable model matricies in =es ce 


& 3.44 are arranged as follows: 


FA MATRIX 
ES 
PS 
Fle O 
Fle - 
Fte 
Fte 
0) Fa 
7 Fa 
Fr 
Jade 
GA MATRIX 
Gst 
Gst 
Gle @) 
Gle 7 
Gte 
Gte 
Ga 
@) Ga 
as Gr 
Gr 
HA MATRIX 
Hst 
Hst 
Hle 
Hte @) 
Hte ~ 
Hle 
Hle 
0 Ha 
~ Ha 
Hr 
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APPENDIX G 


Geb ies misOkKS TRANSFER FUNCTIONS 


fee RATE GYRO TRANSFER FUNCTION 


S 
----- + lige @ 
Y(S) 31.7 
ES) S S 
[----- + 1.0] [----- + 1.0] 
167.8 461.7 


ee ocCELEROMETER TRANSFER FUNCTION 


S 
—---- are oe) 
reals ) 225.6 
E(S) S 2 2(.96) 
(Sasoe ee Smee © 
3952453 8959.3 


meee 2NGLE OF ATTACK SENSOR TRANSFER FUNCTION 


Note: Transfer functions input deagrees/second (or degrees) 
and output degrees/second (or deqrees). The state variable 
model for the transfer functions in Eas. 3.46 & 3.47 are 
developed using the procedures outlined in appendix F. 


ley 


SubroOvea ne: 


Deseriot ion: 


Calling sequence: 


APPENDIX H 


SIMULATION PROGRAM “SUEROULINiGg 


jo i Ra |S) 


Reads in steady state flight conditions, Stee 
and rudder commands, and Control sumnaee 
failure parameters. 


CALL FLITE1 (mach,alt,alpha,nz,ncontymeree 
NStp,anp, est nyt see 


Input arguements: None 


Output arguements: 


mach 
alt 
alpha 
az 

ts 


Neone 


NSe, Nstp 


amp 


ESt ft Aa Sstue 


Steady state mach number 

Steady state altitude in feet 
Steady state AOA in dearees 
Steady state normal acceleration. Normally 1.0. 
Sampling time in seconds 


Control number: l=Longitudinal stick 


2=Lateral stick 
3=Rudder deflection 


Control start time, and control stoprcimes 

Nst and nstp are actually the iteration number: 
Start time (secs. )=nst*ts 
Stop time (secs. )=nsto*ts 


Control amplitude given as deflection in inches 


Right and left stabilator failure parameters 
respectfully: No failure = O 
Failure = 1 


they: 


Pmlearoutine: FLITE2 
Operation: Reads in hasic airframe matricies from the 
Shy emcee 1 le. 
mmeima arguments: CALL FLITE2(Fx,Gx,Hx,Dx,Fyz,Gyz,Hyz,Dyz, 
Nise, No xen Nox Ntyz, Novyz, Nnyz, Ndyz ) 


Miput arguments: None 


Seeput arguments: 


mee, GX, Hx, Dx Longitudinal state matricies 
mz, Gyz, Hyz, Dyz Lateral-directional state matricies 
mex, Nox, Nhx,Ndx Two-dimensional vectors giving the 
tiyz, Noyz,Nhyz,Ndyz number of rows and columns of the 
respective matricies. Fxample, 
Nfx(1) = Number rows in Fx 
Nfx(2) = Number columns in Fx 


roo 


Subroutine: 


Description: 


LONLAT 


Generates the LONG and LATD matricies 
(Eqs. 3.38 and 3.39) hased on tne contre 
surface failure parameters. 


Computes the unimpaired input matrix, GrdO, to 
the modified aiframe equations. The GmO matrix 
will be used in the reconfiguration algorithm. 


Calling Arguments: CALL LONLAT(rstf,lstf,ifail,Gx,GyzZ pene 


LATD,Gm0O,Nlong,Nlatd,Ngm0,Nax,Nayz) 


Input arguments: 


oe Oe) Ao) eee ol 


ifail 


Gx, Gy Zz 


Nox; Noyvz 


Right and left stabilator failure parameters 
respectfully. (See subroutine FLITE1) 


Failure flag: No failure = 0 
Failure = 1 


Longitudinal and lateral-directional input 
matricies. 


Two dimensional vectors giving the number of 
rows and columns of the Gx, and, Gyz 
matricles. 


OUtpPue arguments: 


LONG, LATD 


GmO 


Nlkong,. Niaed 
NgmO 


Matricies described in Section II1I.B. whe 
split the control surface deflections into 
right and left hand sides. The LONG and LATD 
matricies reflect control surface failure or 
damaqde. 


The unimpaired input matrix to the modified 


airframe equations (Fq. 3.40). Gm0O is 
composed using the unimpaired LONG and LATD 
matricies. It 1s subsequently used in the 


reconfiguration algorithm. 


Two dimensional vectors givina the number of 
rows and columns in the respective matricies. 
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Subroutine: AIRDAT 


Description: Computes the air data inputs to the control 
laws hased on standard atmosphere conditions. 


Calling arguments: CALL AIRDAT(mach,alt,temp,rho,psi,a, 


ri) 
Input arguements: 
mach Steady state mach number 
alt Steady state altituce in feet 
Sueput arguements: 
temp Standard atmosphere temperature in degrees 
rankine 
sc lgto, Standard atmosphere density in 


ip-see 6 / f£t**4 


psi Standard atmosphere pressure in lb/ft**2 
a Sonic velocity in ft/sec 

qe Dymamicwpsessure in lbh/ft**2 

iggal Pressure ratio = qc/psi 
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Subroutine: | MOLEo 


Description: 


Composes the mocified airframe matricies (Eas. 


3.40 and 3.41) 


Calling Arguments: 


CALL MODEQO(Fx,Gx,Hx,Dx, PyzZ,GyvzZ) Hy ae 


LONG, LATD, Fm, Gm, Hm, Dm, N£x, Nox, Naxx, Nae 
Nfyz,Ngyz,Nhyz,Ndyz) 


Input arguments: 


nx , Gx , imo. , bx 
BYZ, GY2Z, fyZ,. by 2 


LONG, LATD 


Nex, Nex, Nax,, tex 
Nivyz,.Noyvz, Nhyz, Ney z 
Nilonag;, Niatd, 


Output arguements: 


en, Gin, Him, Dm 


Nim, Ngm, Nhm, Ndm 


Longidudinal and lateral-directional 
basic airframe matricies (Fas: 3a 
- 3.35) Read in from sulrouemee 

Pai Gree: 


3.38 
Generated in subroutine 


LONG and LATE mat rierec Joes 
and 32.39) 
LONLAT 


Two dimensional vectors giving the 
number of rows and columns in the 


.respective matrix. 


Modified basic airframe matricies 
(Fos. 3240 anaes. 41) 


Row and column vectors as described 
as described ahove 
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Subroutine: SENSOR 


Description: Composes the sensor matricies (Faqs. 3.47 and 
3.46) 


Calling arguments: CALL SENSOR(Fs,Gs,Hs,Nfs,Ngs, Nhs) 
Input arguments: None 


Output arguements: 
ms, Gs, Hs Sensor matricies 


Nfs, Ngs, Nhs Two dimensional row and column vectors 


Subroutine: ACTU 


Description: Composes the actuator matricies (Faqs. 3.42 and 
Segers 


Calling arguements: CALL ACTU(Fa,Ga,Ha,Nfa,Nga, Nha) 
Maput arguments: None 


Beeputc arguments: 
Fa, Ga, Ha Actuator matricies 


Nfa, Nga, Nha Two dimensional Row and column vectors 
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Subroutine: PLANT 


Description: 


Composes the basic airframe plus actuator 


matricies (Eqs. 3.44 andes. 


Calling arquments: 


Input arguments: 


Pin. Ging Fm, Dm 


Fa, Ga, Ha 


Nfm, Ngm, Nhm, Ndm 
Nfa, Nga, Nha 
OULDUE arguements: 


Boy Gh ato 
Nip, Nop. Nap 


CALL PLANT(Fm,Gm,Hm,Dm, Fa,Ga,Ha,Fp,Gp, Hp, 
Nfm,Ngm, Nhm,Ndm,Nfa,Nga,Nha,Nfvo,Nap,Nhp) 


Modified basic airframe Matyr vere 
(Eqs. 3.44 and 3.45) Generateamag 
Ssubrouttnew ObEoOe 


Actuator matricies (Eqs. 3742 anawemmees 
Generated in subroutine ACTU. 


Two dimensional row and column ve@teme 


Airframe plus actuator matricies 


Row and column vectors 
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Subroutine: LAWS 


Description: Computes the function gains (App. A) 


Computes the coefficients in the control law 
matricies 


Composes the control law matricies 


Salling arguments: CALL LAWS(alpha,nz,psi,aqc,ri,ts,Ac,Bfc, 
Ba Ce, bite, be. Nac, Note, Noe, Nec, Nafc, Ndc) 


Pmyout arguments: 


alpha Steady state angle of attack in decrees 
mz. Steady state normal acceleration = 1.0 
psi standard atmosphere pressure at steadv state 


Slee Geoimaliny £t* «2 


qc dynamic pressure in 1lb/ft**2 
ri pressure ratio 
ts sampling time in seconds 


Smeput arguments: 


Ac, Bfc, Bec CogenG ee eawematracies “Oras. 3.29 and 3.30) 
me, Dic, De 


Nac, Nbfc, Nbc Row and column vectors 
mee, Ndfc, Ndc 
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Subroutine: VGAIN 


Description: Computes the variable gain matrix based on 
control surface impairment. The SUbroemema 
is disigned to use reconfiguration algorithm 
Introduced ane ket. 


Calling arguments: CALL VGAIN(ifail,ifix,Gm,Gm0,GAIN,Nom, 
NgmO, Ngain) 


1peUe —arouments: 


ifail Control surface failure flag set in subroutine 
-LONLAT; No failure = O 
Failure = l 
sia tao Fix parameter. The GAIN matrix will be computed 
fo the impaired aircraft if ifgx — 1 ane 
ifail = 1. Ifix is read interactively by sub- 


routime 2b Chae 
Gm Airframe control matrix. From subroutine MO@aes 
GmO Unimpaired airframe control matrix. Generated 
in subroutine LONLAT. Note erro noemerc ae 


control surfaces are damaged or failed, Gm=Gro0Q. 


Nom, Ngm0O Row and column vectors 


Output Arguements: 


GAIN Variable gain matrix described in section 
i A Se 2 
Ngain Row and column vector 
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Bioroutine: EXPINT [ Ref.5] 


Description: Computes both the matrix exponential 


Fops*ts 
Aps = e 
and the integral 
ik Fpos*s 
f e ds 


Oo 


Reference 5 gives the method used to perform 
the computations. 


Calling arguements: CALL EXPINT(Fps,Nfps,Aps,Naps,Duml, 
NEM le nes hom, Dime) 
Input arguments: 


Fps Airframe plus actuator plus sensor transfer 
matrix. Generated in ARCRFT 


Nfps Two dimensional row and column vector 
ts Samplina time 
10p Print parameter: O = Do not print results 


Foss ts 
Otherwise print Fps, Ts, e 


Dum2 Vector Or working space for computations 


Output arguments: 


Aps Matrix exponential described above 
Naps Two dimensional row and column vector 
Duml 'The integral of the matrix exponential. To he used 


in computing the discrete control matrix Bps. 


43x10 43x43 43x10 
Bps = Duml x Gps 
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The following subroutines perform the matrix operations 
necessary to compose the A/F-18 system matricies. 


Subroutine: MULT [ Ref .5] 
Description: Performs the matrix multiplications. 
Calling arguments: CALL MULTIA,NA, B/NEB, Ge @ 


Input arguments: 


A, B Matricies to be multiplied 


Na, Nb Two dimensional vector giving the number of rows 


and columns an Ajeanease 


Output arguments: 
C Product of A and B, “€ = G38 


Ne Row and column vector 


Subroutine ADD and SUBT are called just as MULT. 
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Subroutine: NULL [Ref. 5] 

Description: Generates a null matrix 

Calling arguments: CALL NULL(A,NA) 

Maput arguments: 

Na Two dimensional vevtor giving the number of rows 
ane cOLunns so fe eneumt) ie maer) xX, 


Na) = =Nutmper of rows 
Na’ 2Je=oNumboer Of columns 


@arput arguments: 


A Null matrix having dimension Na(1) x Na(2)} 


Subroutine: JUXTR [Ref. 5] 


Description: Juxtaposes two matricies by row 


Calling arguements: CALL JUXTR(A,Na,B,Nb,C,Nc) 


Meaput arguements: 
rh 2B The two matricles to be juxtaposed by row 


Na, Nb Row and column vectors 


Output arguments: 


cc The juxtaposition of A and B C= --- 


Ne Row and column vector 
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Subroutine. VJUXTC ) i Reta on 


Description: Juxtaposes two matricies by column 


Calling arguments: CALL JUXTR(A,Na,B,Nb,C,Nc) 


Input arguements: 
A, B The two matricies to be juxtaposed by column 


Na, Nb Row and column vectors 


Output arguments: 


Il 
Pod 
se 


‘& The juxtaposition of A and B C 


Ne Row and column vector 


Subroutines: ADD and MULT 
The calling sequences for these subroutines are exactly as 


described for those matrix operations above. 


Le) 


Smmoroutine: OUTPUT 


Description: Outputs a matrix with rows and columns 


numbered. 
Calling Arguments: CALL OUTPUT(A,Na(1),Na(2),°A *) 


Powe arguments: 
A Matrix to he output 


hel), Na(2) Mumisen Of roews ana number of columns 
in A respectfully 


A { Matrix name. Must he placed in quotes 
and occupy four spaces. 


Output arguments: 


i7eeA matrix 1S written to the FAL& RESULTS file 


iL Dal 


APPEND 


F/A-18 EXEC PROGRAM 


&TRACE OFF 

CLRSCRN 

GLOBAL TXTLIB VALTLIB VFORTLIB IMSLDP NONIMSL CMSLIB 
GLOBAL LOADLIB VFLODLIB 

FILEDEF 06 TERMINAL 

FILEDEF O01 DISK F/A-18 DATA Al 

FILEDEF 02 DISK F/A-18 RESULTS (LRECL 132 
FILEDEF 03 DISK OPTMATD DATA Al 

FILEDEF 04 DISK OPTPLOT DATA Al 

&TYPE LOADING FA18 SIMULATION 

LOAD FA18 ORACLS1 (START 

& EXIT 


Is 


AXE 


e. 


APPENDIX J 


moccuil CONDI Li@hssAND STABILITY AND CONTROL DERIVATIVES 


FLIGHT CONDITIONS 


MACH = 0.6 
ALPHA = 2.6184 
PCTHLC = 20.496 
DLHTD = 0.68914 
DLAD = 0.0 
Peerl, = 3.0968 
DTEFL = 3.2647 
DRUDL = 0.0 
FLARUD = 0.0 

Se) — 052206 
pervY= 0.1 220E+06 


ALT 


Iho Us One ab 
ve THRE 


DRHD 
DRAD 


DLEFR 
DTEFR 
DRUDR 


DELS 
WAIT 
XIZZ 


T 


B 


ol 


LONGITUDINAL DERIVATIVES 


EOCUO 2c 
Sis) oi ies 
20.496 

0.68914 

0.0 


| | 


ae, 

celbeered 6 (U2 10a 
-64114E-02 
-0.45250E-01 


= 


0 
e) 
0 
6) 


LATERAL-DIRECTIONAL DERIVATIVES 


= —0O,.13257E-0O1 

a) / L265E-0] 
=m), 2090 

Bee) 0.39729E-03 
SB = 0.0 
TH = 0.14257 
S = -0.43034E-03 
LF = -0.412942E-03 
TF = 0.17¢674E-03 
= —0.2437 

De= 0.0 

=m. 27567 

=m .6/478E-02 
A = 0.86683E-03 
pee 0.5291 8E-02 
LF = 0.0 

AT = -0.21027E-02 
Tee— 0.0 


=0 291 25E—01 
O50 

On OO76 

— eee ol 
-0.43808F-02 
Oy ZO56E—0g 

= 0.0 

= 0.39478E-02 

= 0.24948 


(igs 3: 


VTRUE = 646.42 
GAMA = 0.0 
Dee = 21227700 
REX Z = = 1027 20 
MU = -0.12988E-04 
MW = -0.11331E-01 
MO = -0.59346 
MWD = -0.34049F-03 
MDSB = 0.0 
Moi = 0.0 
MDS = -0.24151E-02 
MDLF = -0.5230E-03 
ibe ea O. 32/7 3-0 | 
NV = 0.69481E-02 
NVD = 0.0 
NR = -0.21705 
NP = -0.15254E-01 


NDA = 0.77305E-04 
NDR = -0.46503E-03 
NDLF = 0.0 

NDHT = 0.57361E-04 
NDTF = -0.32655E-02 


ELIGHT CONDITION PARAMETERS 


DEFINITION 


MACH NUMBER 

ALTITUDE, MSL 

TRUE AIRSPEED 

ANGLE CF ATTACK 

TOTAL ENGINE THRUST 

FLIGHT PATH ANGLE 

LEFT THROTTLE POSITION COMMAND (1O04=MAX AB) 
RIGHT THROTTLE POSITION COMMAND(TOOR=MAX AB) 
LEFT HORIZCNTAL STAB. DEFLECTION (+ TE DOWN) 
RIGHT HORIZONTAL STAB. DEFLECTION (+TE DOWN) 
LEFT AILERON DEFLECTION (+TE DOWN) 

RIGHT AILERON DEFLECTION (+TE DOWN) 

LEFT LEADING EDGE FLAP DEFL (+LE DOWN) 

RIGHT LEADING EDGE FLAP DEFL. (+LE DOWN) 

LEFT TRAILING EDGE FLAP DEFL. (+TE DOWN) 
RIGHT TRAILING EDGE FLAP DEFL (+TE DOWN) 

LEFT RUDDER DEFLECTION ( + TE LEFT) 

RIGHT RUDDER DEFLECTION (+TE LEFT) 

FLARED RUDDER ( PA/TO OKLY ) 

SPEEDBRAKE DEFLECTION 

LONG. CENTER OF GRAVITY LOCATION IN PCT MAC 
AIRCRAFT WEIGHT 

MOMENT Cr INERTIA ABOUT X-BOOY AXIS 

MOMENT OF INERTIA ABOUT Y-BODY AXIS 

MOMENT OF INERTIA ABOUT 2-BOOY AXIS 

XZ PLANE PRODUCT Cf INERTIA 
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RAD 
PERCENT 
PERCENT 
DEG 

DEG 

DEG 

DEG 

DEG 

DrG 

DEG 

DEG 

DEG 

DEG 

DEG 

DEG 
PERCENT 
LB 
SLUG=Fi2 
SLUGS 
SLUGSE tz 
SLUG =F 


FLIGHT CONDITIONS AND SAMPLING TIME 


AIRCRAFT 


Q0.2618E+01 


E+O5 ALPHA 


aC 


+01 TS 


ORM 


AND FAILURE PARAMETERS 


CONTROL PARAMETERS, 


AMPLITUDE 
SAB 


SLART Sire 
RIGHT 


Sereno tx K 


F/A-18 RESULTS FILE 
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